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ABSTRACT 


The  sputtering  process  has  been  investigated  by  simulating 
the  sputtering  of  single-crystal  copper  with  1-7  keV  argon.  A 
digital  computer  was  used  to  build  the  crystal,  bombard  it,  and 
move  crystal  atoms.  Four  mechanisms  were  observed  which 
cause  surface  atoms  to  sputter.  An  atom  is  sputtered  when 
(1)  it  is  squeezed  out  of  the  surface,  (2)  it  is  scooped  out 
when  another  atom  strikes  its  inner  hemisphere,  (3)  it  Is  ejected 
when  an  atom  passes  behind  it,  and  (4)  it  is  knocked  out  by  a 
second  layer  atom  which  is  moving  outward.  Nearly  all  sputtered 
atoms  were  surface  atoms.  Second  and  third  layer  atoms  were 

t  ' 

sputtered  only  for  ion  energies  greater  than  5  keV.  They  'were 
sputtered  by  mechanisms  similar  to  the  surface  atom  mechanisms. 
"Silsbee  chains"  were  observed  to  be  directed  into  the' Crystal, 
and  momentum  focising  was  observed  to  cause  sputtering  only 
when  It  occurred  in  close  packed,  surface  rows.  Outward 
directed  chains  were  not  observed.  Sputtering  deposit  patterns, 

w  ....  . .  .  ....  ,  .  .  v. 
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sputtering  ratios,  and  Sputtered  atom  energy  distributions  were 
obtained  for  ( 10Q>,  {110} ,  and  (111)  surfaces.  All  data  , 
compared  favorably  with  experimental  data. 
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1.  INTRODUCTION. 

Material  is  lost  from  a  substance  undergoing  ionic  bombard¬ 
ment.  This  is  sputtering.  Many  authors  have  tried  to  define 
sputtering  using  terms  such  as  "(atoms) ..  .ejected  or  knocked 
out",  "emission",  "ionic  erosion",  "disintegration",  and  "break¬ 
down."  A  reader  who  blithely  accepts  one  of  these  definitions 
in  the  course  of  his  reading  may  conclude:  This  definition 
describes  the  mechanism  of  sputtering.  This  is  an  illogical,  If 
not  erroneous,  conclusion;  the  definition  of  sputtering  must  be 
derived  from  the  mechanism,  not  the  mechanism  from  the 
definition.  The  sputtering  process  has  been  investigated 
utilizing  one  of  two  general  mechanisms,  ejection  or  emission. 
Ejection  implies  an  immediate  or  direct  reaction  release  of  an 
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atom  from  a  crystal;  emission,  a  delayed  or  indirect  reaction 

process  which  results  in  the  release  of  an  atom. 

( 1) 

In  1923,  Kingdon  and  Langmuir  7  bombarded  thoriated 
tungsten  with  various  ions  in  a  glow  discharge  tube.  This  was 
a  special  case  of  sputtering  since  the  thin  surface  film  of 
thorium  on  a  tungsten  substrate  was  sputtered  rather  than  the 
tungsten  Itself.  The  results  of  this  experiment,  that  the 
sputtering  ratio,  atoms  removed  per  incident  ion,  increased 
with  increased  ion  mass  and  Increased  ion  energy,  qualitatively 
suggested  an  ejection  mechanism.  A  few  years  later,  Von 
Hippel  and  Blechschmidt^  proposed  a  theory  which  described 
sputtering  as  an  evaporation  of  surface  atoms,  an  emission 
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mechanism.  Earlier*  Von  Hipped 3)  had  found  by  spectroscopic 
means  that  at  least  some  sputtered  atoms  were  in  an  excited 
state.  The  sputtering  theory  showed  that  atoms  In  the  region 
of  impact  could  rapidly  acquire  thermal  energy  if  the  kinetic 
energy  of  incident  ions  was  converted  to  thermal  energy  at  the 
target  surface.  If  an  atom  acquired  a  sufficient  amount  of 
thermal  energy*  It  would  then  evaporate  from  the  surface,  some 
atoms  evaporating  while  in  an  excited  state.  This  theory  was 
improved  by  Townes ^  in  1944,  who  calculated  a  sputtering 
flux  and  the  number  of  atoms  evaporated  per  incident  ion. 

Key w ell ^  in  1954,  used  neutron  diffusion  theory  to  approxi¬ 

mate  atomic  interactions  within  the  crystal,  a  new  approach  to 
the  theoretical  investigation  of  sputtering.  Direct  application 
9,f  statistical  methods  to  sputtering  was  made  by  Harrison^ ) 
who  envisioned  the  Interaction  of  two  distribution  functions 
(the  crystal  lattice  and  the  ion  beam).  These  models*  based  on 
statistical  methods,  implicitly  accept  ejection  type  mechanisms. 

One  of  the  most  important  contributions  to  the  study  of 
sputtering  was  made  by  Wehner^  in  1953  •  In  the  first 
sentence  of  his  paper*  Wehner  stated*  "The  most  widely 
accepted  sputtering  theory  is  the  evaporation  theory...",  but 
the  deposit  patterns  of  single  crystal  sputtering  which  he 
obtained  showed  pronounced*  high  density  areas  or  "spots".. 

There  was  now  strong  evidence  for  a  momentum  transfer  , 
process,  and  further  development  of  the  evaporation  theory 
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ceased.  Shortly  after  Wehner’s  findings  were  reported, 
Henschke^  proposed  a  theory  of  sputtering  based  solely  on 

t 

classical  collision  theory,  treating  normal  and  oblique  incidence 

•* 

sputtering  separately.  Oblique  incidence  sputtering  could  be 
explained  by  two  body  collision  processes.  The  incident  ion 

penetrated  the  surface  layer  with  little  or  no  interaction,  was 

« 

relected  outward  by  second  layer  atoms,  and  ejected  a  Surface 
atom  by  striking  it  on  its  inside  hemisphere*  His  concept  erf 
normal  incidence  sputtering  required  many-body  collisions  In 
which  the  ion  waa  eventually  reflected  outward  to  sputter  v.  • 
surface  atoms*  This  theory  was  plausible  for  oblique  Incidence 

ft 

sputtering,  however,  the  case  for  normal  incidence  sputtering 
requited  that  the  ion  be  reflected  inside  the  crystal,  {For  an 
ion  more  massive  than  the  target  this  requirement  cannot  be 
met.) 

Bllsbee  noted  that  in  any  discussion  of  momentum  trans¬ 
fer  effects  in  sputtering,  the  geometry  of  the  crystal  should 
be  consideredi  successive  collisions  in  a  crystalline  structure 
might  be  influenced  by  the  structure  itself.  His  calculations 
showed:  in  a  two-body  collision,  an  atom  has  a  departure  angle 
with  respect  to  an  axis  specified  by  a  close  packed  direction. 
In  subsequent  collisions,  the  departure  angle  will  decrease 
If  the  energy  is  low  enough,  and  the  collision  sequence  occurs 
along  a  close  packed  row  of  atoms.  This  concept  of  "momentum 
focusing",  "BUsbee  chains",  or  "focusbfts"  was  considered  by 


many  investigators ,  excepting  Wahair^^.,  to  be  a  fully  satis¬ 
factory  explanation  of  the  spots  in  deposit  patterns#,  But  the 
deposit  pattern  is  only  one  of  many  observable  features  of 
single  crystal  sputtering. 

The  sputtering  ratio  is  another  Important  characteristic 
of  the  sputtering  process.  If  a  theory  or  model  Is  to  explain 
sputtering,  both  patterns  and  sputtering  ratios  must  be 
explained  qualitatively  and  quantitatively.  Almen  and  Bruce, 
(AB)  measured  sputtering  ratios  of  a  variety  of  metals 

using  N,  Ne,  Ar,  Kr,  and  Xe  ions  over  a  5-65  keV  range  of 
bombardment  energies.  (Although  poly  crystalline  specimens 
were  sputtered  ,  the  qualitative  results  are  the  same  as  those 
obtained  for  single  crystal  sputtering.)  They  noted  that  the 
sputtering  ratio  was  an  increasing  function  of  bombardment 
energy  provided  the  ion  mass  was  greater  than  the  target  maos 
Sputtering  ratios  showing  a  maximum,  or  of  k  slowly  varying 
nature,  were  observed  in  all  cases  in  which  the  ion  was  lighten 
Evidence  such  as  this  coincides  with  any  gross  conception  of  an 


ejection  mechanism. 

i 

Single  crystal  copper  was  sputtered  by  argon  at  inter¬ 
mediate  energies  by  Magnuson  and  Carlson  (MC)^i2^  and  South¬ 
ern,  Willis,  and  Robinson  (SWR)*13*.  MC  measured  sputtering 
.ratios  ( using  1-10  ksV  argon)  for  the  (ill),  (100)  and  (110) 
surfaces  and  found  that  sputtering  ratios  decreased,  in  this  ' 
order  of  surfaces.  SWR  sputtered  single  crystal  copper  with 
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1-5  keV  argon  and  measured  sputtering  ratios,  but  of  greater 
significance  was  the  quality  of  their  deposit  patterns  which 
clearly  showed  the  presence  of  only  certain  spots s 

i  ‘ 

(111)  surface  3-(110),  l-(lll)  but  not  seen  due  to 

beam  aperture* 

(110)  sir  face  l-(HO),  2-(100) 

(100)  surface  4-(ll0),  1-(100) 

The  work  of  AB,  MC,  and  SWR  Indicated  that  the  sputter¬ 
ing  ratio  is  at  least  a  function  of  bombardment  energy  (or 
momentum)  and  of  mass  ratios.  However,  consideration  of 
spot  patterns  indicates  that  the  sputtering  process  Itself  Is 
also  highly  dependent  on  the  gross  and/or  surface  geometry  of 
the  crystal  and  on  Ion  penetration  depth. 

Ion  penetration  of  solids  has  been  Investigated  experiment¬ 
ally  and  by  computer  simulation.  Piercy,  McCargo,  Brown, 
and  Davies (Ik)  investigating  channeling  of  various  heavy  ions 
in  monocrystalline  aluminum,  found  that  the  ion  penetration 
distance  increased  in  the  order  (111),  (100),  and  (110)  beam 
orientations.  This  is  the  same  order  as  the  sputtering  ratio 
decrease  found  by  MC.  A  number  of  studies  of  ion  penetra¬ 
tion  have  also  been  made  using  computer  simulations  (Oen, 
Holmes,  and  Robinson^  *5)  f  Robinson  and  Oen^1^,  Harrison, 
Deeds,  and  Gay^1?)). 

The  most  apparent  correlation  between  the  results  of  ion 
penetration  and  sputtering  studies  is  the  concept  of  trans¬ 
parency.  defined  by  Flult,  Rok  and  Kistemaker^8^.  Both 
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experimental  and  simulation  investigations  of  ion  penetration 
indicate  that  channels  exist  in  certain  crystallographic  directions 
of  a  crystal.  Results  of  normal  incidence  sputtering  experi¬ 
ments  indicates  that  sputtering  ratios  are  lower  for  surfaces 

i  •  *  ..  ...  .  .  • '  .  *  t ,  '  ■  } 

in  which  these  channels  are  parallel  to  the  beam.  One  might 
initially  suspect  that  sputtering  ratios  will  be  low  if  ions  and 
target  atoms  are  confined  to  these  channels  j  momentum  is 
directed  into  the  crystal  rather  than  laterally,  reversal  of 
momentum  to  produce  sputtering  not  occurlng.  Accordingly, 
some  incident  ions  should  see  a  transparent  surface  rather 
than  a  wall  of  atoms. 

The  obvious  complexities  of  the  sputtering  process  may  be 
investigated  separately  by  digital  computer  simulation.  This 
of  course  involves  the  selection  of  a  gross  mechanism  (ejection 
or  emission)  and  a  model  (hardsphere  or  otherwise).  If 
factors  such  as  the  interatomic  potential  function  and  its 
parameters  are  known,  then  the  use  of  the  right  model  in  the 

computer  may  be  expected  to  simulate  the  actual  sputtering 

* 

process.  Explanation  of  deposit  pattern  characteristics  and 
sputtering  ratios  will  necessarily  follow  from  an  exact 
simulation. 


2.  FURPOSE  OF  THE  INVESTIGATION, 


The  belief  that  the  formation  of,  spots  in  sputtering  deposit 

•*  ’  .•  il  ; *  S  ft  * 

patterns  is  solely  a  consequence  of  momentum  focusing  was 
widely  accepted  for  a  number  of  years.  The  Importance  o? 
momentum  focusing,  not  only  in  spot  formation,  but  in  the 
entire  sputtering  process  is  now  believed  to  be  much  less  than 
originally  thought  ( 18—25) .  This  investigation  was  undertaken  to 
explore  the  single  crystal  sputtering  process  using  a  digital 
computer  to  simulate  a  copper  crystal  which  is  being  bombarded 
with  a  beam  of  argon  ions.  Primary  effort  was  directed 
towards: 

a.  Finding  the  mechanisms  which  result  in  atoms  being 
deposited  in  the  spots  or  high  density  regions  of  the 
deposit  pattern. 

b.  Determining  the  correlation  between  crystal  structure 
and  sputtering  ratios  and  patterns . 

c.  Investigating  the  surface  binding  en ergy  and  its  effeet 
on  deposit  patterns. 
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3 .  MODEL. 

A.  A  space  lattice  is  established  whose  sites  represent 

the  equilibrium  positions  of  copper  atoms  In  a  face  centered 

cubic  crystal.  Atomic  spacing  is  that  determined  by  x-ray 

o 

crystallographic  studies  (for  copper,  a  -  3.615  A),  The 

"  ■  i 

potential  function,  with  parameters  for  copper-copper  inter¬ 
action,  is  the  Born-Mayer  type  Gibson  Number  Two^^,  ;  The 
interaction  between  argon  and  copper  is  similarly  described 
using  a  Born-Mayer  potential,  V{  r)  -  eA+Br,  but  the  numerical 
values  of  the  parameters  A  and  B ,  for  this  function  are  those 
determined  by  Harrison,  Carlston,  and  Magnuson 
from  a  study  of  secondary  electron  emission.  The  range  of  both 
potential  functions  is  eroded  at  rQ,  one  half  the  nearest  neigh¬ 
bor  distance  such  that  the  potential  and  force  go  to  sero  for 
atoms  whose  centers  are  farther  apart  than  2rQ.  Use  of  this 
eroded  form  of  the  potential  functions  permits  establishment 
of  the  crystal  in  a  stable,  static  state  without  use  of  an 
attractive  potential.  This  equilibrium  stats  is  destroyed  only 
when  an  atom  moves  from  Its  lattice  site. 

The  lattice  as  a  whole  is  restricted  by  the  current  computer 
program  to  perfect  form.  Defects  such  as  interstitial  atoms 
cannot  be  used,  however,  trregulartties  in  the  surface  layer  are 
introduced  by  using  two  types  of  variation  from  a  perfect 
surface  *  vacancy  and  stub  conditions.  The  perfect  or  regular 
surface  is  identical  to  a  full  lattice  plane.  In  the  vacancy 
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configuration  the  atom  which  would  normally  be  hit  first  by 
the  ion  is  removed  from  the  perfect  surface.  The  stub 
eurface  has  an  atom  added  on  the  perfect  surface  in  a  stable 
position  and  adjacent  to  the  target  atom.  Additional  atoms 
can  also  be  added  on  the  surface  or  removed  from  the  perfect 

i 

surface  to  provide  a  variety  of  random  irregularities* 

B.  Each  ion  of  the  beam  is  approximated  by  a  single, 
neutral  argon  atom  whose  velocity  vector  intersects  the  plane 
of  surface  atoms  at  an  Impact  point.  (The  term  ton  Is  used 
throughout  this  thesis  for  the  Incident  particle  to  avoid  con¬ 
fusion  with  target  atoms.) 

Successive  runs  are  made  using  different  Impact 
points  to  simulate  use  of  the  entire  beam  area.  The  use  of 

this  simple  model  for  the  ion  beam  is  subject  to  four 

•  > 

conditions  t 

1 .  The  results  obtained  using  one  random  set  of  Impact 
points  are  essentially  the  same  as  those  obtained 
using  any  other  random  set  of  impact  points. 

2.  The  majority  of  beam  ions  are  neutralised  before 
reaching  some  arbitrary  boundary  which  defines  the 
surface. 

3.  Prior  to  neutralisation,  the  path  of  a  beam  ton  Is 
not  appreciably  affected  by  the  surface  potential. 

h.  A  crystal  region  is  in  an  equilibrium  state  each  time  • 
a  neutralised  ton  Impacts  In  that  region. 


2 1 


t 

i 


Ea  ch  of  these  conditions  has  been  satisfied  either 

'  f .  *rt  i  r  ■;  ••  \ 

in  the  model  itself  or  by  calculation.  A  full  discussion  is 
contained  in  Appendix  A. 
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4.  DYNAMICS 


The  force  functions  for  copper-copper  and  argon-copper 
interactions  are  derived  from  the  respective  potential  functions. 
Atom  position  and  velocity  are  then  determined  using  Newton’s 
Second  Law  of  Motion.  Normally,  one  would  determine  these 
values  by  numerical  integration  methods  over  small  time  inter¬ 
vals,  however,  this  process  uses  too  much  computer  time. 
Special  methods  approximating  integration  have  been  developed 
which  are  used  to  minimize  computer  run  time  yet  maintain  a 
good  approximation  to  the  integration. 

A.  Qay (2®)  developed  an  iterative  method  of  solution  which 
is  similar  to  that  used  by  Gibson,  Go  land,  Milgram,  and 
Vineyard^  26 ) .  Both  Gay’s  and  GIbson*a  methods  use  the 
technique  of  replacing  differential  quantities  with  finite 
differences  to  approximate  the  equation  of  motion,  F/m  -  dv/dt 
>/  -t  by  F/m  «  Av/  At.  Uaing  this  equation  and  the  relation¬ 
ship#  Ax  »  v&t  and  v  ■  (v  ♦  vQ)/2,  equations*  (3-1)  and  (3-2) 
are  obtained. 

Av  ■  y  -  vQ  -  [f  UQ)/ra]  At  (3-1) 

(3-2) 

The  known  values  F(x  ) ,  vc,  x0,  and  At  in  these  equetloae 
are  usually  associated  with  a  time  t  -  tQj  x  and  f,  the 
unknown  values,  are  to  be  evaluated  using  these  equations  at 
time  t  •  T(t0,  At)* 

Gibson's  technique  associates  the  valve  of  x  with  a  time 

o 


Ax  -  x  -  x. 


-  ^[F  At/2  ♦  ^At 


t  m  t  .  but  v  is  associated  with  a  time  t  ■»  t  -  At/2, 
oxo  o  ov  .  '  o 

Current  values  of  x  and  rare  alternately  computed  at  Intervals 
At  using  equations  (3-1)  and  (3-2). 

vvtQ  +  At/2)  -  v0  (tQ  -  At/2)  +  ^F(xQ  ,/m]At  . 
x(tQ  .At)  =  x{tQ)  +  |[F(x0)/mjAt/2  +  v(tQ  +  kt/2)^ 

The  advantage  of  this  technique  over  that  of  evaluating  both 


variables  at  time  t  -  tQ  +  A  t  is  that  a  continuous  smoothing 
of  the  values  of  x  and  v  occurs  with  each  successive  computa¬ 


tion. 

Gay  recognized  the  need  for  some  sort  of  smoothing 
procedure  but  considered  that  an  averaging  process  for  equation 
parameters  rather  than  computed  results  would  result  in  a 

t 

better  approximation.  His  method  replaces  the  evaluated 
force  function  F(xQ)  in  these  two  equations  with  an  arbitrary 
force  function  f ( x)  which  is  linear  in  some  Interval  to  be 
determined.  Consider  the  following  two  step  cycles 

Move  an  atom  from  xQ  to  x^  using  computed  values  of 
F(x0)  and  vQ  to  solve  equations  (3-1)  and  (3-2).  Compute 
F(xj)  and  average  this  value  with  F(xQ).  Now,  move  the 
atom  from  x^  to  xg  using  the  averaged  force,  F(xc,  x^) 
and  v0  to  solve  the  equations  using  the  same  constant  value 
of  At.  Let  f  now  be  defined  in  the  interval  (x0,  Xg),  such 
that  f(xQ)  -  F(x0,  x^) .  Equations  (3-1)  and  (3-2)  may 
now  be  written: 


v(t  +  At)  ■  v(t)  +  Jf(xo)/mJ  At 
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xtt  +  At)  -  x{t)  +( 

The  values  of  v  e»d  x  are 
which  f  approximates  IT  to  some  arbitrary  degree,  therefore, 
in  every  complete  cycle  or  timeetep  a  new  interval  and  new 
function  f  are  defined* 

A  determination  of  which  of  these  two  methods  Is  better 
must  be  based  on  some  standard  or  condition.  In  many  cases 
the  accuracy  of  numerical  solutions  suck  as  these  may  be 
compared  to  the  exact  value  one  obtains  by  integration.  But, 
there  are  no  truly  exact  results  in  a  simulation  since  a  physical 
model  is  in  the  computer  rather  than  the  computer  being  used 
as  a  means  of  solving  equations.  However,  if  the  model  Is 
physical,  physical  laws  must  be  satisfied,  and  the  energy  gain 
or  loss  due  to  the  mathematical  approximations  must  be  small. 
This  condition  of  energy  balance  has  been  used  by  Harrlso©^*^ 
who  found  that,  in  simulations  involving  atomic  interactions  such 
as  channeling  or  sputtering,  use  of  Gay’s  method  resulted  in  a 
better  energy  balance. 

B.  The  value  of  these  approximation  techniques  may  be 
substantially  reduced  by  injudicious  selection  of  a  numerical 
value  for  At.  Too  large  a  value  invalidates  the  approximation, 
but  too  small  a  value  increases  computer  running  time.  The 
program  used  in  this  simulation  of  sputtering  incorporates  not 
only  Gay’s  approximation  method  but  also  a  procedure  which  he 
developed  for  automatic  adjustment  of  At.  Its  value  Is 


jjf(x0)/m]At/2  +  v(t)j  A  t  ' 
determinable  In  any  interval  in 
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calculated  using  the  relationship  ft  «  fr/v  where  At  is  replaced 
here  by  ft  to  indicate  that  its  value  varies  with  each  ti meat ep . 

m.  vat  - 

The  velocity,  v  is  the  magnitude  of  the  velocity  ,pf  tlje  most 

'  ‘  *  .«»■  *  '  ■*.  -+• 

energetic  atom  in  that  timestep,  and  ^r  is  used  here  aa  an 
input  parameter,  the  times tep  multiplier .  This  parameter  is 
assigned  the  computer  variable  name  DTI  and  is  defined*  the 
maximum  displacement  of  any  atom  In  any  timestep  shall_.be  the 
numerical  value  of  DTI  in  appropriate  units  of  length.  Gay 
considered  the  value  of  &t,  consequently  that  of  DTI,  to  be  a 
function  of  ion  energy,  impact  point,  and  start  point,  but  , 
he  was  unable  to  find  the  relationship.  Trial  and  error  methods 
were  used  to  find  optimum  values,  and  reasonably  good  results 
were  obtained  using  these  values.  Johnson  ^ 29)  reported that 
the  energy  balance  was  maintained  to  within  $%  for  all  his  runs* 
The  difficulties  in  choosing  a  value  for  DTI  have  been  , 
largely  eliminated  in  this  study.  Some  problem  areas  such  as 
start  point  and  impact  point  dependence  were  avoided  by  always 
positioning  the  ion  tangent  to  the  first  atom  It  would  hit. 

This,  at  worst,  will  result  In  a  constant  errqr  for  which 
adjustments  may  be  made  (this  has  not  been  necessary).  The 
positioning  procedure  is  described  in  Appendix  5.  The  depend- 
ence  on  ion  energy  is  Inherent  to  the  DTI  computation  prpcqss 
and  Is  easily  explained  once  the  clutter  of  impact  point  a,p4 
start  point  dependence  has  been  cleared  away.  Use,  of  thf 
DTI  process  excludes  all  variable  dependence  except  the  inverse 
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proportionality  to  velocity,  but  the  interval  inwhioh  the  fores 
function  la  considered  linear  la  exactly  the  same  widt&  as  the 


numerical  value  of  DTI*  Accordingly,  a  straightforward 
analysis  la  made  to  show  the  interdependence* 

The  force  function,  F(r),  is  expanded  about  an  arbitrary 
point  r*  using  a  Taylor's  series  for  a  function  of  one  variables 
F(r)«F(r*)  +^F  I  (r-r*)  *  1  VV I  (r-r»)2+***l  (r«*#) 


F(r)»F(r*)  +^F 

Tr"  r* 


2!*r2 


Since  F(r)  is  derived  from  the  potential  function  V(r)'  •  ^A+Br, 

°  ’  li  F(r)  is  necessarily  of  exponential  form  arid ‘  m  'B^F 

“id  the  series  is  rewritten  in  the  following  form?"  ‘ 

F(r>  ^F(r*)  [l  +  B(r-r*)j|  +  ^F(r*)  Bn  (r-r*l°|| 

The  quantity  (r-r*)  is  now  defined  as  a  variable  displacement 
J^r,  and  equation  (3-3)  may  be  written?  F(r)  -  f  (  S  r)  +  g(  5  P) 

where  t[$  r)  and  g(  J  r)  are  the  terms  in  braces  in  the  force 

( 

function  expansion.  The  function  f(<fr)  Is  assumed  to  be  the 
linear  function  used  in  Gay's  approximation  method}  the  function 
f  (  ^ r)  is  to  approximate  F(r),  Thus,  it  is  required  that  the 
ratiotfjg(  ^r)/f(  &  r)  be  some  fraction  less  than  one* 

Let  an  assumption  be  made  that  the  ratio  is  much  less 
than  one,  but  that  the  actual  value  is  the  fractional  deviation' 
from  linearity,  a  variable  whose  value  is  to  be  specified.  Then, 
since  B  is  a  negative  constant  (from  the  potential,  function), 
the  truncation  error  for  the  series  is  easily  found,  and  the 


quantity  Bdr  Is  known  for  every  specified  value  of  fractional 

1  1  •  '  '  t  V  ?  O'j  ’{.I'lijESnoct'tOqOfXq 

deviation.  But  B  is  constant,  and  once  the  fractional  devia- 

r  ri  1  ■  -  -1  i.if  •  yi 

tion  is  specified,  <Tg  is  fixed.  Consequently,  specified  to  the 

.  r :‘y  "•  -v  -•  .  •  .  ;  '"‘iT:;.-  *  ,  i  .  J  5 ls.‘.-0 1 ^ S mir ;** 

fractional  deviation  from  linearity,  there  is  a  function  f(r) 

■  '■  ■'  •••'  '  '  •'  ’  =»t  f  .?.£  3;^.  - 

•which  is  a  valid  approximation  to  the  function  F(r)  in  every 

■  1  •  . :  -  ■'  ,-t r: -I  - •  -  ^.a:r 

.interval  (r,  r+5r)  where  St  is  fixed.  The  quantity  Sn  is 

thus  the  value  to  be  used  for  DTI,  and  £t  **-.£r/v  -  DTI/v. 

;  ‘  [  ■  \  2^  ■  i  ■’"••:!  .  >  -  H 

It  is  now  assurred  that  at  will  always  he  the  maximum 

possible  value  which  minimizes  computer  running  time  without 

‘  ’  -"V:  •  ■  •  -  \  '  'i  ’'O.uB 

Introducing  errors  due  to  the  non-linearity  of  the  true  force 

_ ...  :•  -  7  - *  { J 

function.  One  would  expect  therefore,  that  the  energy 

’  ■  .  '  •  1  ■  i. 

balance  would  be  maintained  to  the  same  order  as  the  linearity 

"  :  ■  :  "  \ . ~  •  ■  ""  '  '  '■  V'  '  "  ••  ’  •• T)  1 

of  the  force  function.  Further  consideration  will  show  that 
the  energy  balance  is  maintained  to  a  considerably  better 
degree  than  that  anticipated.  First,  since  DTI  is  the  maximum 
displacement  of  the  most  energetic  atdnu  essentially  all  other 
atoms  will  be  displaced  a  distance  less  than  DTI*  Thus* ‘  the 
interval  in  which  the  forte  function  f  (r }  Is  used. ih  smaller* 
and  f  Is  an  even  better  approximation  to  EV  Second,  for  a 
given  fractional  deviation  from  linearity,  the  value  is 

total  displacement,  not  2,  j,  or  a  component  displacement  j 
•  and  the  velocity  used  in  the  determination  &t  2  DTI/v  for 
each  tlmeetep  ie  the  total  velocity  y,  not  the  cfompohient 
'  Velocity  vjV  The  result  of  these  considerations  may.:  be 
Summarised  by  a  calculation  for  the  maximum  displacement  of 
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an  atom  along  the  component  of  its  displacement  IfeCtCr* 

£Xj  -  Vj&t,  but  Vj  £  *  4  ▼maxand  4t  -  DTI/vmax.  It  follow  e 
that  dxi  -  ^  (DTi;/vmax)  C  DTI. 

These  analyses  appear  to  provide  a  tidy  solution  to  the 
time  optimization  and  energy  balance  problem.  This  would  be 
true  except  for  the  complications  introduced  by  the  eroded 
form  of  the  potential  function.  Consider  a  situation  in  which 
one  atom  is  separated  from  another  by  a  distance  (2rQ  +  d) 
where  d  is  some  distance  less  than  the  total  distance  O  that 
the  atom  will  move  in  the  next  times tep.  When  the  atom  Is 
moved  the  distance  D,  as  shown  in  figure  lt  the  force  on  tike 
atom  at  its  new  position  is  the  same  force  it  would  have  tf 
it  had  moved  only  (D-d)  units.  Since  the  force  is  aero  for 
r>2rc,  the  model  assumes  that  the  atom  had  a  velocity  v  * 
(D-d)/6t  (which  Is  smaller  than  Its  actual  velocity  v  »  D/&  t.  , 
and  a  smaller  kinetic  energy  is  computed  using  this  smaller 
velocity.  This  difficulty  was  originally  foreseen  by  Gay  who 
Included  corrections  to  the  force  calculations  for  such  situa¬ 
tions.  Both  Levy ^30)  and  Johnson^ 2^)  improved  the  methods 
of  corrections,  and  recent  Improvements  in  the  force  calcula¬ 
tions  have  been  made  by  Harrison^  D , 

The  continuous  improvement  of  the  dynamics  section  of 
the  program  has  resulted  in  a  model  in  which  not  only  has  the 
computer  runhing  time  been  minimised,  but  of  greatest  signifi¬ 
cance,  for  nearly  all  of  tike  700-800  combinations  of  surfaces, 
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Reports  of  sputtering  studies  often  include  sd  hoc  formula¬ 
tions  of  mechanisms  which  are  used  to  explain  deposit  patterns, 
sputtering  ratios,  and  sputtered  atom  energy  distributions* 

The  sputtering  process  Is  discussed  in  this  thesis  In  terms  of 
mechanisms  which  have  been  observed  In  »•>*«  rimi^atlesi  to  mm attar 
atoms*  The  observance  of  these  mechanisms  la  an  advantage 
peculiar  to  a  simulation.  Bach  crystal  atom  (and  the  Ion)  must 
be  Identified  by  number,  at  least  for  purposes  of  computer 
calculations.  Since  every  atom  is  identified,  Its  complete  track 
can  be  plotted  and  labeled.  The  tracks  of  selected  atomic 
then  be  superimposed  to  show  a  complete  set  of  interactions . 

The  mechanisms  which  were  observed  to  cause  an  atom  to 
sputter  are  considered  the  prime  observable  quantity  of  the 
simulation. 

The  results  of  the  simulation  are  discussed  for  each  crystal 
surface.  The  determination  that  an  atom  haring  energy  £  is 
sputtered,  was  made  using  a  probablllty-of -sputtering  function, 
P(B)  »  1  -  e“®/®b  ethane  1%  is  the '.assumed  .binding,  energy* 

The  value  of  P  Is  compared  with  a  random  number  R  having  a 
value  between  0  and  1*  If  PVR,  the  atom  Is  considered 
sputtered.  (This  method  of  selection  has  been  eomiortd  with 
one  using  a?s€ep  function  in  which  all  atoms  with  B>£^  are 
assumed  sputtered.  A  selection  method  ustef  aether  grobabtt- 
ity-of-aputterfng  function,  P(i)  »  i  -  hss  also  been 
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tested.  No  differences  were  observed  in  the ^puttering 
.  patterns  based  on  the  atom  selection  method*)  .  ,yf,_ 

Mechanisms  observed  to  cause  sputtering  are  discussed,, 

_•  * 

first  for  each  surface.  Deposit  patterns ,  sputtering  ratios, 
and  sputtered  atom  energy  distributions  are  then  discussed*: 
The  deposit  patterns  shown  are  superposition*..  Qf „  those 
obtained  by  sputtering  the  crystal  separately  with  regular , , 
vacancy,  and  stub  surface  conditions.  This  practice  follows, 
that  used  by  .Johnson^  y?ho  found  that  pattern  features  - 
were  difficult  to  recognize  when  viewing  separate  patterns 
from  each  surface  condition.  The  energy  distributions  are 
also  superpositions  of  the  three  surface  conditions*  The  use 
of  superposition  is  considered  to  most  realistically  simulate 
the  condition  of  the  crystal  surface  at  various  times  during 
sputtering.  There  are  strong  arguments  however,  for  using 
only  the  regular  surface  condition.  Sputtering  ratios  are 
not  averaged  (which  would  correspond  to  the  superposition 
of  other  data)  since  a  probability  factor  for  the  existence 
of  each  surf  ace  condition  would  need  to  be  determined*  In¬ 
stead*  the  simulation  sputtering  ratios  for  vacancy  and  stub 
surface  conditions  are  included  and  discussed  briefly  only  for 


completeness.  (The  (110)  surface  was  sputtered  only  at  1 

* 

and  3  keV  using  all  three  surface  conditions. ) 

Crystallographic  nomenclature  ie  used  generloally  in  the  r 
discussion  of  spots*-  Reference  to  gn  OdcU  spot  d@ee.j9ot 


Imply  {hkl^ ejection.  The^hld^dlrectione^are  with  respect  to 
a  right  handed  coordinate  system  in  which  the  crystal  ta 
described*  The  (hkl)  surface  is  contained  in  the  x-s  planet 
the  2  direction  is  into  the  dryetai.  The  three  crystals  used 
in  the  simulation  are  shown  in  figures  2-4 1  all  contain  15d 
atoms.  This  crystal  sise  has  been  found  sufficient  to  contain 
nearly  all  energetic  collisions  for  ion  energies  up  to  7  -*V. 

.  a  .  __ 

A.  The  (111)  surface. 

The  (111)  surface  was  sputtered  normally  with  argon  at 
1-5 1  7(  10,  20,  and  40  keV.  (Although  the  potential  function 
is  considered  valid  only  for  ion  energlee  less  than  about  7keV, 
the  higher  energy  runs  were  used  to  search  for  additional 
mechanisms.)  Four  mechanisms  were  found.  These  Were 
especially  evident  at  the  lower  ( 1-5  keV)  ion  energies.  Three 
of  these  are  classified  surface  mechanisms  since  only  surface 
atoms  were  found  to  participate  in  the  sputtering  event.  The 
fourth  mechanisms  Is  a  deep  mechanism  In  which  an  inward 
moving  atom  Is  reflected  and  then  initiates  a  sputtering  event. 
In  all  cases,  the  dominant  mechanisms  were  found  to  be 


Except  for  head-on  collisions,  the  Ion  le  scattered  by  the 
target  (atom  2)  with  a  component  of  momentum  parallel  to 
the  surface  as  well  as  a  much  larger  perpMdtohtar  component 
Into  the  crystal.  The  effect  of  this  combination  of  momenta 
le  to  drive  a  nearest  neighbor  (a.  n. ),  for  #xampl*  atom  6  hi 
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figure  5,  into  an  apparent  (ll2>  channel  formed  bjrStome  7, 

8 ,  and  9  •  The  channel  is  apparent  since  it  terminates 
abruptly  thus  causing  atom  6  to  drive  atoms  7  and  8  Into  the 
crystal  and  atom  9  outward.  Since  atom  6  passes  nearly 
directly  behind  atom  9,  the  impulse  is  more  normal  to  the 
surface  than  parallel  to  it.  This  is  termed  a  mole  mechanism 
since  one  atom  burrows  between  two  layers  parallel  to  the 
surface  to  sputter  an  atom  in  the  outermost  layer.  Atoms 
sputtered  by  this  mechanism  are  almost  always  a  next  nearest 
neighbor  (n.n.n.)  to  the  target  and  located  in  or  near  the 
sextant  defined  by  the  impact  area.  They  usually  sputter 
with  greater  than  10  ev.  The  target  always  receives  the 
majority  of  energy  transferred  by  the  ion,  but  it  is  normally 
driven  into  the  (111)  trigonal  array  of  atoms  directly  behind 
It.  These  three  atoms  act  as  buffers  and  dissipate  the 
target's  energy  and  momentum  into  the  crystal.  The  target 
is  reflected,  but  it  does  not  retain  sufficient  energy  to 
sputter  a  surface  atom  nor  sputter  itself. 

The  second  and  third  surface  mechanisms  occur  when  the 
impact  parameter  is  about  one  third  of  an  atom  radius  or 
less.  In  both  mechanisms  a  n.n.,  atom  6,  is  struck  by  the 
ion  such  that  the  parallel  and  perpendicular  momentum  compon¬ 
ent  magnitudes  which  it  acquires  are  of  the  same  order.  If 
the  perpendicular  component  is  greater,  atom  8  scoops  a  n.n.n., 
atom  9«  relatively  high  up  In  the  surface  plane  of  atoms  and 
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into  their  outer  hemispheres.  These  atoms  are  densely  packed 
in  the  surface  and  act  as  a  rigid  reflector  for  atom  9*  It  is 
found  that  atoms  sputtered  by  the  scoop  mechanism  hare 
ejection  angles  nearer  the  normal  than  would  be  found  ill  the 
absence  of  reflection  in  the  surface.  If  the  perpendicular 
momentum  component  of  atom  6  is  nearly  the  same  or  less 
than  the  parallel  component,  it  squeezes  atom  9  against  Its 
surface  neighbors.  The  squeezing  causes  the  surface  plane  to 
warp  and  atom  9,  10,  or  11  is  sputtered.  Atoms  sputtered 
by  either  the  scoop  or  squeeze  mechanism  usually  sputter  with 
less  than  10  ev.  Atoms  sputtered  by  combinations  of  the 
three  surface  mechanisms  are  found  to  sputter  with  higher 

I 

energies  than  when  sputtered  by  a  single  mechanism. 

The  fourth  mechanism  requires  that  an  atom  be  reflected. 
The  atom  which  is  most  often  apt  to  be  reflected  Is  the 
target.  If  the  impact  parameter  is  greater  than  one  half 
the  atom  radius,  the  target  is  driven  towards  an  edge  of  the 
trigonal  array  behind  it  rather  than  into  It.  It  wlU  penetrate 
this  edge  if  it  has  sufficient  energy,  but  tt  will  lose  meet  of  the 
energy  during  the  penetration.  The  reflection  occurs  from 
third  layer  atoms  and  the  target  atom's  reversed  momentum  Is 
transferred  to  a  second  layer  atom.  The  second  layer  atom 
will  then  sputter  one  of  two  surface  atoms  at  energies  upwards 
of  5  ev.  At  lower  Ion  energies  the  target  atom  cannot 
penetrate  the  array  edgs  and  is  reflected  from  the  second 
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layer j  it  may  or  may  not  sputter  a  surface  atom,  but  it 
rarely  sputters  itself. 

When  the  crystal  was  sputtered  at  7  keV,  the  ion  was 

t 

found  to  penetrate  well  into  the  second  layer.  This  is  a 
significant  increase  in  ion  penetration  depth  since  5  keV  ions 
are  found  to  penetrate  only  to  the  order  of  half  an  atom  radius 
The  effect  of  deeper  ion  penetration  is  to  cause  a  lateral 
compression  of  the  second  layer  with  accompanying  warping. 
Second  layer  atoms  which  are  squeezed  outward  will  sputter  a 
surface  atom  but  are  rarely  found  to  sputter  themselves.  At 
the  higher  bombardment  energies  tested  (10,  20,  40  keV),  the 
number  of  second  layer  atoms  which  sputtered  increased  nearly 
proportionally  to  the  increase  in  ion  energy.  The  sputtering 
process  at  these  higher  energies  appears  to  be  predominantly 
by  the  three  surface  mechanisms  described  for  lower  energy 
sputtering.  Mechanisms  peculiar  to  high  energy  sputtering  have 
not  been  observed  although  this  may  be  due  to  the  limited 
energy  range  for  model  validity. 

All  (111)  simulation  deposit  patters  showed  the  character¬ 
istic  features  of  three  (110)  spots  and  a  central  (111)  spot. 
The  2  keV  pattern  definitely  showed  the  presence  of  three 
additional  spots  which  were  found  in  the  regions  in  which 
streaking  occurred  at  other  energies.  These  patterns  are 
shown  in  figures  6-21.  The  appearance  of  a  hexagonal  pattern 
at  2  keV  was  a  surprising  result,  but  it  Is  r.ot  a  unique 


occurrence  In  sputtering  studies .  Anderson  and  WeJmer^ 
found  a  hexagonal  pattern  for  (111)  copper  sputtered  by 
mercury.  The  expected  trigonal  pattern  was  found  for 
energies  up  to  400  ev,  a  hexagonal  pattern  appeared  at  400 
ev,  and  it  disappeared  as  the  ion  energy  increased.  Very 
recently,  Robinson  and  Southern^^)  have  found  additional 
spots  near  positions  for  (111)  gold  sputtered  with 

4  keV  argon. 

The  pattern  features  vary  slowly  with  binding  energy 
provided  its  value  is  restricted  to  the  range  2. 50-3 •  59  ev. 

The  upper  limit  is  the  sublimation  energy  for  copper  (33)j 
above  it  the  patterns  appear  to  deteriorate.  This  was  most 
evident  at  a  3  keV  bombardment  energy.  A  set  of  patterns 

'j 

at  this  bombardment  energy  is  shown  in  figures  8-14  for  a 
binding  energy  range  of  1 . 50-5*50  ev.  The  pattern  deteriora¬ 
tion  Is  most  evident  in  the  region  of  the  (111)  central  spot. 
Figures  11  and  12  show  the  transition  from  a  well  defined 
central  spot  at  3«00  ev  to  complete  deterioration  of  the 
central  spot  at  4*00  ev. 

The  half -intensity  width  of  the  (110)  spots  is  estimated 
from  pattern  data  and  from  numerical  data  to  be  about  11 
degrees.  The  width  of  the  central  spot  is  an  unreliable 
datum  and  was  not  determined. 

Sputtered  atoms  were  always  found  to  be  surface  layer 
atoms  for  sputtering  at  5  keV  or  less.  The  atoms  which 
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were  sputtered  most  frequently  appeared  to  naturally  group  by 
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energies  Into  two  categories ;  atoms  with  perpendicular  energy 

.  -  .  *  '  •  .  '  •  • ' '  '  r;  j  j I ' O  i 

greater  than  3  but  less  than  10  ev  and  atoms  with  perpendicular 
energy  greater  than  10  ev.  The  higher  energy  atoms  were 

•  -  f  *  "'V'  ;.j" 

found  to  be  equally  distributed  between  the  three  (110)  spots 

and  the  central  spot.  They  were  also  found  in  the  narrow 

sectors  defining  the  streaks  and  at  distances  from  pattern 

center  corresponding  to  the  distances  of  the  (110)  spots. 

This  distance  was  found  to  be  0.7  units,  the  distance  measured 

for  the  spots  in  the  patterns  of  SWR.  The  lower  energy 

atoms  were  usually  found  in  the  central  spot  which  explained 

its  sensitivity  to  binding  energy. 

The  correlation  of  sputtering  frequency  and  crystal 

location  of  sputtered  atoms  was  found  using  f requency-location 

diagrams.  These  are  shown  in  figures  22-2^  for  regular 

•  * 

surface  sputtering.  Similar  diagrams  have  been  made  for 
sputtering  the  vacancy  and  stub  surfaces.  It  was  found  that 
a  1  keV  ion  does  not  cause  frequent  sputtering  of  n.n.n’-'s  or 
n.n.'s  in  the  sextant  defined  by  the  impact  area.  However, 
at  2  and  3  keV,  n. n.n.'s  are  found  to  be  sputtered  most 
frequently,  and  at  h.gher  energies,  both  n. n.n.'s  and  n.n.'s 
are  frequently  sputtered.  The  three-fold  relationship  between 

I 

the  frequency  of  sputtering,  location  of  the  atom  with  respect 
to  the  target,  and  the  ion  energy  may  be  clarified  by  consider¬ 
ing  the  sputtering  mechanisms. 
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At  1  keV  and  an  impart  parameter  of  about  three- 1 berths 
an  atom  radius,  the  ion  is  found  to  penetrate  less  than  half 
an  atom  radius.  At  energies  up  to  5  keV  the  penetration  is 
only  slightly  greater.  It  was  found  that,  although  the  target 
is  always  driven  into  the  crystal,  the  n.n.  which  initiates  Idle 
majority  of  sputtering  may  not  be  driven  far  enough  into  the 
crystal  to  provide  the  scoop  mechanism.  If  it  is  not  driven 
inward  a  sufficient  distance  to  scoop  the  n.n.n. ,  its  energy 
will  be  propagated  in  the  surface  along  close  packed  rows 
which  originate  at  the  n.n.n.  This  is  not  a  rcfocuson.*1  At  an 
early  point  in  the  propagation  the  surface  will  have  warped 
sufficiently  to  cause  an  atom  to  sputter.  The  sputtering  is 
a  result  of  a  squeeze  mechanism  but  not  one  directly  involving 
the  n.n.  As  the  ion  energy  is  Increased  the  scoop  and  mole 
mechanisms  are  more  apt  to  occur  although  the  squeeze 
mechanism  is  found  to  occur  at  all  ion  energies.  The  transi¬ 
tions  from  dominance  of  one  mechanism  to  another  with 
increasing  ion  energies  may  be  Inferred  from  the  profiles  shown 
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In  figure  31  for  atoms  which  sputter  most  frequently. 

A  relatively  constant  sputtering  frequency  is  observed  for 
atoms  102  and  132,  but  the  percentage  of  high  and  low  energy 
sputa  varies  for  atom  102  which  is  a  n.n.n.  The  percentage 
does  not  vary  for  atom  132  which  is  neither  a  n.n.n.  nor  a 
n.n.  Atom  102  sputters  most  frequently  at  low  ion  energies 
by  the  squeeze  mechanisms,  at  hl^ier  ion  energies  by  the  scoop 
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and  mole  mechanisms.  These  relationships  have, been  confirmed 
by  the  atom  track  displays.  They  also  have  shown  that  atom 

■  *'  V  •  r>‘'  "■  .*  ^  i  O:.-  .  ".VibtS'i  rftCXi  J?  l  r.JO 

132  is  sputtered  by  atom  102  through  a  squeeze  mechanism  ' 

*  ■  ■  -  ■  .  ■  -V\  r  u- 

alone  at  low  ion  energies  but  in  conjunction  with  tiie  mole 

mechanism  at  higher  ion  energies.  The  profile  shown  for  atom 

*• 

86  is  the  result  of  simple  reflection  from  second  layer  atoms. 

% 

Atoms  87  and  101,  which  are  not  in  the  impact  sextant  but 
are  sputtered  frequently,  are  sputtered  by  the  squeeze 
mechanism . 

Figure  32  shows  profiles  of  atoms  sputtered  by  the  deep 
mechanism.  At  low  ion  energies  the  target  reflects  at  an 
Oblique  angle  from  second  layer  atoms  rather  than  penetrating 
the  layer,  and  it  enters  one  of  the  apparent  <(ll2^  channels. 
Atoms  25  and/or  55  are  then  sputtered  directly  by  the  target. 

:  t  . «  i  ■ 

(Although  this  is  similar  to  a  mole  mechanism  It  is  a  distinct 
mechanism  since  the  target  must  be  reflected  in  this  case  but 
not  in  the  former  case.)  At  higher  Ion  energies  the  target 
will  penetrate  the  second  layer  and  Is  reflected  from  third 
layer  atoms.  It  does  not  channel,  but  it  causes  a  second 
layer  atom  to  sputter  surface  atoms.  Accordingly,  atom  55  is 
sputtered  by  the  deep  mechanism  more  often  than  atom  25  Is 
sputtered  by  a  channeled  target.  Atom  40  or  70  is 
sputtered  by  a  squeeze  mechanism,  usually  in  conjunction  with 
the  sputtering  of  atoms  25  and  55* 

These  results  for  the  sputtering  of  a  regular  surface 


have  been  found  to  be  generally  applicable  to  sputtering  Vacancy 
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and  stub  surfaces.  Equal  numbers  of  atom*  from  the  three 
surfaces  are  found  in  the  spots  and  in  the  streak  region* • 

It  is  also  found  that  these  atoms  are  visually  the  same  from 
all  three  surface  conditions .  They  dif f er  only  in  energy  since 
it  is  the  ion  penetration  distance  which  indirectly  determines 
the  magnitude  and  direction  of  momentum  transferred  to 
surface  atoms.  The  Vacancy  surface  sees  an  ion  which 
penetrates  deeper;  this  is  analogous  to  an  ion  of  greater  energy* 
In  the  stub  condition,  the  ion  does  not  penetrate  the  surface; 
it  transfers  energy  to  the  stub  atom.  The  stub  then  assumes 
the  role  of  the  ion  impinging  on  the  surface ,  but  the  stub 
appears  as  an  ion  with  less  energy. 

Figure  33  shows  the  sputtering  ratios  for  each  surface 
condition  as  a  function  of  ion  energy;  assumed  binding  energy 
is  a  parameter.  A  statistical  variation  of  +  20%  of  the 
sputtering  ratio  has  been  assumed.  This  would  normally  be  an 
unreasonably  large  deviation  for  laboratory  results ,  but  it  is 
considered  conservative  for  the  numerical  results  from  a 
simulation  in  which  many  parameters  are  unknown.  The  actual 
value  of  sputtering  ratios  obtained  for  the  regular  surface  is 
generally  low  by  2-3  atoms/ion  for  a  binding  energy  of  3.5  er. 
Further  investigation  showed  that,  in  a  few  cases,  surface 
atoms  3  planes  below  the  Impact  point  would  be  sputtered  with 
about  4-5  ev.  The  average  Increase  in  the  (111)  sputtering 
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ratio  over  the  1-7  keV  range  Is  estimated  to  be  no  greater 
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than  2.0  atoms/ion  for  a  binding  energy  range  of  1.5  to  3.5 
ev.  This  range  is  shown  by  the  dotted  lines  in  figure  32a. 
Vacancy  and  stub  sputtering  ratios  show  erratic  behavior  oyer 
the  energy  range  studied. 

Energy  distributions  of  sputtered  atoms  are  shown  in 
figures  35-39*  The  preponderance  of  atoms  with  energies  less 
than  5  ev  is  caused  by  the  lack  of  an  intrinsic  surface  binding 
energy.  The  sputtering  selection  process  considers  all  atoms 
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having  E  greater  than  1  ev.  Thus  the  P(0.50)  energy,  using 
an  assumed  binding  energy  of  3.5  ev,  Is  2.1*2  ev.  This  is 
especially  evident  for  the  1  and  5  keV  distributions.  If  this 
region  is  disregarded,  the  peaks  which  occur  at  5-8  ev  for 
all  ion  energies  are  assumed  to  be  the  maximum  for  each 
distribution.  Secondary  maxima  appear  at  16  ev  for  1  keV 
sputtering,  16  ev  for  2  keV,  and  at  11*  ev  for  3  keV  sputter¬ 
ing.  At  higher  ion  energies,  a  secondary  maximum  may  be 
present  at  1*7  ev  for  5  keV  sputtering.  The  number  of  atoms 
sputtered  in  the  simulation  is  too  small  to  make  a  definite 
statement  regarding  maxima.  Differences  are  usually  measured 
by  one,  at  most,  two  atoms. 

5.  The  (100)  Surface. 

The  (100)  surface  was  sputtered  normally  with  1,  3,  5, 
and  7  keV  argon.  Mechanisms  observed  to  cause  sputtering 
were  found  to  be  nearly  Identical  to  the  surface  mechanisms 
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dtscusaed  for  (111)  sputtering.  A  mechanism  directly  ' 
comparable  to  the  deep  mechaniem  was  not  observed.  Instead, 
it  was  found  that  the  scoop  and  mole  mechanism#  are  enhanced 
by  the  presence  of  the  (110)  channels  parallel  to  the  surface. 
Surface  atoms  which  are  driven  into  these  channels  do  not 
have  to  burrow  between  the  first  and  second  layer  In  order 
to  sputter  nearby  surface  atoms.  Additionally,  It  was  found 
that  these  atoms  will  frequently  sputter  themselves  by 
reflection  from  second  layer  atoms.  The  aqueese  mechanism 
was  observed  to  be  the  most  effective  sputtering  mechanism 
at  3  and  5  keV.  At  1  and  7  keV  ion  energies,  variations  of  the 
scoop  and  mole  mechanisms  were  dominant.  The  low-high  and 
middle  ion  energy  dependence  observed  for  the  dominant 
mechanism  is  directly  related  to  the  shape  of  the  impulse  by 
which  the  tion  transfers  energy  to  a  n.  n. 

At  low  Ion  energies,  the  impulse  is  sufficiently  broad  that 
the  n.n.  is  directed  into  the  edge  of  the  square  of  atoms 
behind  it.  The  n.n.  will  enter  the  (110)  channel  without 
immediately  sputtering  a  surface  atom.  When  the  impulse  Is 
narrow,  the  peak  force  is  not  necessarily  greater.  The  target 
has  received  most  of  the  energy  given  up  by  the  ion,  and  the 
ion-n.n.  Impact  parameter  Is  larger.  The  n.n.  will  then  scoop 
or  &  surface  atom  causing  it  to  sputter.  At  high  ion 

energies  (7  keV)  the  ion  was  found  to  penetrate  at  least  Into 
the  second  layer  where  It  was  deflected  towards  but  not  into 


43 


i 


a  (110)  channel.  This  was  observed  to  cause  secot 


rer  atoms 


to  be  scooped  up  and  sputter  surface  atoms.  Second  layer  atoms 
are  also  squeezed  by  their  neighbors  and  may  sputter  through 
the  vacancy  left  by  a  sputtered  surface  atom. 

Sputtering  at  3  and  5  keV  was  found  tp  be  caused  by  the 
spattering  of  surface  atoms  along  the  surface  as  the  ion  was 
being  reflected  from  second  layer  atoms.  The  squeeze 
mechanism  was  observed  to  sequentially  sputter  atoms  along 
a  close  packed  surface  row.  The  occurrence  of  a  mole 


mechanism  was  conspicuously  rare{  even  at  5  keV  the  energy 

i 

propagation  was  clearly  restricted  to  propagation  parallel  to 

i 

I  •* 

the  surface  layer  of  atoms  and  to  propagation  into  the  crystal. 

■  .  *1  —  .  •  <_  ■ 

The  only  definite  momentum  reversals  observed  were  for  the 


target  atom  and  the  ion.  «  j 

I 

The  characteristic  features  of  (100)  sputtering  deposit  1 

.  »’ 

patterns  are  four  (110)  spots  and  a  central  (100)  spot.  The 
patterns  are  usually  outlined  by  a  hypo  cycloid-shaped  haze)  ( 

the  (110)  spots  form  the  cusps  of  the  hypocyclold.  Simula¬ 
tion  patterns  were  found  to  show  these  characteristic 
features  quite  well  for  3  and  3  keV  sputtering.  The  3  keV 
pattern  (figure  41)  Is  a  very  good  likeness  to  one  at  2.5  keV 
reported  by  SWR.  The  (110)  spot  distances  from  pattern 
center  were  found  to  be  1.0  TC  unit  which  corresponds  to 

e 

the  distances  in  the  2.5  keV  experimental  pattern.  It  was 
further  determined  from  numerical  data  that  the  (110)  spots 
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w«r«  actually  elliptical  rather  titan  circular.  This  latter 
characteristic  is  undoubtedly  the  result  of  distortion  inherent 
In  flat  plate  collection  of  sputtered  atoms  ^he  3  keV 

pattern  was  found  to  remain  essentially  unchanged  as  the 
binding  energy  was  varied  from  2.50  to  3*50  ev.  The  5  keV 
pattern  (figure  42)  was  almost  identical  to  the  one  at  3  fceV 
only  when  a  binding  energy  near  2.00  ev  was  used*  As  the 
binding  energy  was  increased,  the  (110)  spots  became  lees 
well  defined  (figure  43)  and  spots  corresponding  to  { 210} 
became  the  most  prominent  feature.  The  intensification  of 
(210}  regions  was  accompanied  by  a  loss  in  definition  of  the 
hypo  cycloid  outline.  Similar  results  were  obtained  when  spatter* 
ing  at  7  keV.  The  hypocycloid  outline  in  the  pattern  could 
only  be  observed  for  a  binding  energy  near  1.50  ev. 

The  deposit  pattern  for  1  keV  sputtering  (figure  40) 
showed  very  little  similarity  to  those  obtained  from  higher 
energy  sputtering.  Four  spot#  corresponding  to  (21l}>  *»*re 
observed;  four  (liO)  spots  msy  also  bs  defined  but  they  are 
sxtremely  diffuse.  (The  ( 110)  spots  are  best  sssn  in  ths 
point  plot  In  figurs  40s.)  Ths  finding  of  (211)  spots  st 

1  ksV,  and  ths  appearance  of  (210)  spots  for  5  and  7  keV 

•  .... 

Ion  energies  suggested  that  the  pattern  might  be  rotated  45 
degree!  for  low  ion  energy  sputtering.  The  possibility  that  a 
rotation  occur#  ha*  **en  investigated  but  one  has  never 

been  experimentally  observed . 


The  different  (100)  pattern  features  which  appeared  with 
varying  Ion  energies  and  binding  energies  were  a  marked  contrast 
to  the  relative  constancy  observed  for  (ill)  patterns.  (The 
use  of  a  binding  energy  of  3*00  ev  resulted  in  comparable  (111) 
patterns  for  all  ion  energies.)  The  lack  of  constancy  for  the 
pattern  simulations  may  be  explained  as  a  result  of  the  sub¬ 
traction  of  a  binding  energy  from  the  perpendicular  component 
of  a  sputtered  atom's  energy.  (This  Is  done  to  simulate  the 
energy  lost  in  overcoming  the  surface  potential.)  An  atom 
sputtered  with  small  perpendicular  energy,  less  than  15  ev  for 
example,  may  suffer  an  apparently  small  change  in  its  perpendic¬ 
ular  velocity,  but  the  direction  of  its  velocity  may  be  consider¬ 
ably  altered.  The  average  perpendicular  energy  of  atoms  found 
In  spot  regions  is  about  15-20  ev  from  sputtering  the  (100) 
face.  It  Is  about  20-25  ev  for  (111)  sputtering  for  which 
variations  of  up  to  1  ev  binding  energy  have  had  little  effect  on 
pattern  features.  This  oxplt  ition  does  not,  however,  explain 
why  the  pattern  at  3  keV  is  valid  for  a  range  of  binding 
energies.  A  consideration  of  the  atoms  sputtered  may  clarify 
but  not  completely  explain  this  situation.  The  disslmi  arlties 

seen  in  the  (100)  patterns  may  be  directly  Inferred  from  the 

* 

frequency-location  diagrams  for  (100)  sputtering  (figures  46-49). 

When  the  energy  propagation  is  restricted  to  the  surface, 
as  it  is  for  1  keV  sputtering,  the  (111)  surface  dissipate 
energy  through  four  atoms  In  the  momentum  forward  semi- 


circle.  The  (100)  has  only  three  atoms  comparably  located 
to  dissipate  the  energy.  Track  displays  has*  shown  that  for 
1  keV  sputtering,,  atoms  26,  116,  and  132  are  sputtered  by  a 
combination  of  mole  and  squeeze  mechanisms.  The  mole  la 
predominant  in  sputtering  the  atoms  nearest  the  target,  the 
squeeze  mechanism  becoming  predominant  for  surface  atoms 
which  are  located  further  from  the  target.  Atom  26  sputters 
less  frequently  than  atom  116  since  it  is  on  the  opposite  aide  6? 
of  the  crystal  with  respect  to  the  impact  area.  Both  of 
these  atoms  are  found  in  (211)  spots,  but  atom  116  Is  also 
found  in  (110)  regions  of  the  pattern.  Atom  132  la  found  In 
(211)  spots  and  the  (100)  spot{  It  is  sputtered  only  by  the 
squeeze  mechanisms.  When  the  bombardment  energy  Is  3  keV, 
the  deposit  pattern  has  its  expected  characteristics.  Correla¬ 
tion  of  high  atom  density  regions  with  t requency-locatton  data 
showed  that  the  (110)  spots  are  the  result  of  sputtering  of 
n.n.n.'e.  Concurrently,  the  hypocycloid  outline  was  found  to 
be  formed  byt  (1)  the  sputtering  of  the  same  atoms  which 
formed  (211)  spots  at  1  ksV,  and  (2)  the  sputtering  of  this 
n.n.n.'e  In  the  surface.  These  atoms  wars  not  densely 
deposited  in  the  (211)  regions  but  formed  diffuse  elBpeet  whose 
sernt-mlnor  axee  were  along  quadrant  bisectors.  The  appearance 
of  well  defined  (110)  spots  at  3  keV  la  strongly  suggested  by 
the  3  keV  frequency-location  diagram,  figure  47.  Atoms  70 
and  72  were  observed  to  be  sputtered  with  about  the  same 


frequency  as  atoms  55;  56  and  86.  Since  atoms  ?G  and  72 
are  along  a  (100)  axis,  one  might  expect  that  they  would  be 
found  in  (110)  spots,  and  this  has  been  observed.  Atoms  55 
and  56  are  in  positions  relative  to  the  target  such  that  they 
would  sputter  within  adjacent  45  degree  sectors  In  quadrants 
II  and  III.  These  atoms  are  found  in  the  hypocycloid  outline. 

The  absence  of  the  hypocycloid  at  1  keV  and  its  presence 
at  3,  5  and  7  keV  is  readily  apparent  when  one  considers  the 
sputtering  profiles  in  figure  50.  At  1  keV  ion  energy,  n.n.’s 
are  never  sputtered,  however,  they  are  sputtered  with 
increasing  frequet‘;cy  as  the  ion  energy  increases.  Atom  86, 
the  n.n.  in  the  quadrant  containing  the  impact  area,  exhibits 
this  behavior  of  increasing  sputtering  frequency.  Atom  72 
shows  a  sputtering  profile  which  indicates  the  (110)  spot 
formation  at  ion  energies  of  3  keV  or  greater.  The  constant 
nature  of  the  profile  for  atom  102  at  3  keV  and  higher  ion 
energies  is  indicative  of  the  hypocycloid  outline  rather  than 
<2U>  spot  formation.  Atoms  26,  116  and  132  show  a  generally 
decreasing  frequency  of  sputtering  with  high  ion  energy.  This 
is  consistent  with  the  loss  of  (211)  spots  at  higher  ion  energies. 

The  track  patterns  observed  for  5  and  7  keV  (100) 
sputtering  showed  that  atoms  ejected  by  a  squeeze  mechanism 
are  often  ejected  in  directions  opposite  to  those  one  would 
expect  solely  on  the  basis  of  the  location  of  the  atom  with 
respect  to  the  impact  point.  These  .occurrences  were  caused 
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by  two  factorst  (l)  the  atom  was  squeezed  against  Its 
neighbor  and  reflected  outward  and  with  its  parallel  component 
of  momentum  reversed  rather  than  being  strictly  ejected  out¬ 
ward,  (2)  Sputtering  by  neaj  simultaneous  squeeze  and  mole 
mechanisms  results  in  ejection  of  an  atom  where  the  ejection 
direction  is  dependent  only  on  the  impulse  delivered  by  the 
channeling  atom  in  the  mole  mechanism. 

Sputtering  ratios  for  the  ((100)}  sarfaefeaare  dhoWwlB'fifture  51 
The  correspondence  between  simulation  values  and  experimental 
values  is  quite  good  for  the  regular  surface  for  a  binding 
energy  of  3*50  ev.  The  ratios  for  vacancy  and  stub  surfaces 

*v 

show  a  closer  correlation  to  the  regular  surface  for  (100) 
sputtering  then  for  (ill)  sputtering. 

The  energy  distributions  of  sputtered  atoms  are  shown 
in  figures  52-55*  No  specific  maxima  other  than  that  near 
5-7  ev  appeared  to  be  present. 

C.  The  (110)  Surface. 

This  was  the  last  of  the  three  face  centered  cubic  crystal 
surfaces  to  be  sputtered  In  this  simulation.  It  is  purposely 
the  last  to  be  discussed.  First,  experimental  deposit  patterns 
from  (110)  sputtering  show  only  a  large  central  oval  area  for 
bombardment  energies  greater  than  a  few  hundred  ev  (10). 
Second,  a  unique  feature  of  this  surface  is  the  (110)  channel) 
no  other  (hkl)  surface  hi  a  face  centered  cubic  crystal  shows 
^hkl^  channels. ;  The  sputtering  mechanisms  occurring  for  t‘,  3, 


5  and  7  keV  ion  energies  were  moat  readily  observed  for.  . 
sputtering  of  this  surface,  and  they  are  interrelated  with  the 
presence  of  (110}  channels.  The  mechanisms  are,  again,. 
Identical  in  concept  to  those  previously  discussed  for  (111) 
and  (100)  sputtering.  Before  discussing  these  mechanisms 
separately  in  terms  of  individual  crystal  atoms,  the  features 
of  the  frequency-location  diagrams  (figures  56-59)  are 
summarized. 

Sputtering  of  the  (110)  surface  differs  significantly  from 
the  sputtering  of  (100)  and  (ill)  surfaces;  (1)  the  target 
atom  was  the  atom  sputtered  most  frequently  at  all  ion 
energies.  (This  is  in  agreement  with  results  obtained  by 
Levy  ^'.)  (2)  the  number  of  sputtered  atoms  with  high 

energies  were  generally  those  in  the  close  packed  row  contain¬ 
ing  the  target  atom,  but  with  the  target  atom  located  at 
the  center  of  the  row  rather  than  at  the  origin  of  the  row. 
These  findings  are  wholly  consistent  with  expected  results 
when  one  considers  the  sputtering  mechanisms  with  respect  to 
this  particular  surface. 

The  frequent  sputtering  of  the  target  is  made  possible  by 
the  nature  of  the  squeeze  mechanism  but  the  mechanism  does 
not  Itself  cause  the  sputtering.  When  tne  target  is  driven 
into  the  crystal  it  squeezes  the  n.n.  (in  Its  povt)  and  ejects 
It.  Once  the  n.n.  position  is  vacant,  the  target  vfUL  travel 
through  a  large,  potential-free  area  before  striking  a  aeCpnd; 
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layer  atom.  The  target  strikes  this  second  layer  atom  and  is 
free  to  reflect  outward  without  finding  a  surface  atom  directly 
in  its  path.  The  target  then  transfers  most  of  its  parallel 
momentum  to  a  n.n.n.  in  the  row.  The  parallel  impulse  is 
then  propagated  down  the  surface  row.  Warping  of  this  one 
row  causes  additional  atoms  to  sputter.  This  sputtering 
sequence  is  also  observed  as  a  result  of  the  Ion  striking  the 
n.n.  which  is  in  the  quadrant  defined  by  the  impact  area.  A 
rather  surprising  continuation  of  this  sputtering  mechanism 
was  observed  to  occur  in  atom  rows  located  both  above  and 
below  the  horizontal  row  containing  the  target.  Second  layer 
atoms  which  receive  energy  from  the  target  or  ion  are  driven 
Into  (111)  apparent  channels.  They  pass  behind  another  second 

layer  atom  cawsing  it  to  be  ejected  outward  and  strike  two 

*«* 

adjacent  surface  atoms.  Sputtering  ir ,  again,  initiated  In  a 

surface  row.  Sputtering  profiles  for  the  target  (atom  72) 

and  atoms  42,  71,  73  and  102  are  shown  in  figure  60.  Atoms 

42  and  102  are  the  n.n.'s  to  the  target}  atoms  71  and  73  are 

those  which  initiate  sputtering  in  the  horizontal  rows  containing 

■% 

them.  Second  layer  atoms  are  sputtered  more  frequently  as 
the  Ion  energy  is  increased. 

Simulation  depoait  patterns  from  (110)  sputtering  (figures 
61-64)  tend  to  show  more  of  an  oval  outline  than  a  uniformly 
dense  central  oval  region.  The  central  region  can  be  mads 
mors  dense  by  Including  atoms  sputtered  with  less  than  1  sv 


perpendicular  energy  in  the  prob&biiity-of -sputtering  selection 
process.  This  did  not  seem  justifiable  for  (110)  pattern  ... 
production  since  these  small  energy  atoms  were  excludedfrom 
(100)  and  (111)  patterns.  The  1  keV  point  plot  (figure  6la) 
showed  indications  of  (100)  spots  but  these  cannot  be  clearly 
seen  in  the  deposit  pattern  (figure  6lb).  Their  presence  was 
substantiated  by  numerical  data.  The  average  energy  deposit 
data  showed  that  the  average  energy  of  atoms  in  ^100*^  regions 
was  nearly  double  that  of  atoms  found  in  the  central  oval. 

This  is  a  constant  characteristic  of  non-central  spot  regions. 
At  higher  ion  energies  there  was  no  direct  indication  of  t)ie 
presence  of  these  spots  either  by  numerical  data  or  energy 
deposit  data.  A  few  atoms  having  high  energy  were  found  in 
these  regions  but  the  ,  rea  density  was  too  small  to  form  a 
spot  in  the  pattern.  The  set  of  atoms  found  to  form  the 
central  oval  for  sputtering  at  all  ion  energies  was  not  a  well 
defined  group  such  as  has  been  found  for  (100)  and  (111) 
sputtering.  The  atoms  forming  the  (100)  spots  from  1  keV 
sputtering  are  those  which  are  either  n.n.'s  to  the  target 
or  n.n.'s  to  the  atom  which  initiates  sputtering  in  each 
horizontal  row. 

The  sputtering  rati^  for  (110)  sputtering  is  shown  in 
figure,  65.  Agreement  with  experimental  data  from  the 
numerical  and  curve  shape  aspect  was  poorer  for  this  surface 
than  the  other  two  surfaces.  This  is  apparently  caused  by 


52 


{ 110  >  sensitivity  of  the  sputtering  ratio  to  the  A  pare  meter 
In  the  argon— copper  potential  function.  Devy  obtained 

better  curve  shapes  using  A  *=  11.435  whereas  A  »  12. 56 
has  been  used  In  this  study.  This  sensitivity  has  not  bean 
observed  for  {100)  and  (111)  sputtering. 

Energy  distributions  are  shown  in  figures  66-69-  The 
notable  feature  of  (110)  sputtering  is  that  the  energy  of 
sputtered  atoms  is  generally  higher  for  (110)  sputtering  than 
(100)  or  (111)  sputtering. 

D.  Results  c f  Ancillary  Studies. 

During  the  course  of  the  simulation,  it  became  apparent 
that  certain  extensions  of  this  study  should  be  made.  It  was 
particularly  desirable  to  probe  five  areas: 

( 1 )  Sputtering  mechanisms  peculiar  to  high  bombard¬ 
ment  energies  (up  tc  40  keV). 

(2)  Sputtering  of  the  (110)  surface  by  high  energy 
ions  which  are  directed  towards  ( 110)  channels . 

(3)  Sputtering  mechanisms  for  an  ion  heavier  than 
the  target. 

(4)  Sputtering  of  the  (0001)  basal  surface  of  a 
hexagonal  close  packed  crystal. 

(5)  Sputtering  the  (100)  surface  at  ov  ion  energies. 

Investigation  of  the  first  two  of  these  five  areas  were 

limited  somewhat  by  the  restrictions  for  the  copper-copper 
and  argon-copper  potential  function.  They  are  usually  assumed 


valid  up  to  about  10  k eV,  however,  if  small  impact  parameters 
are  assumed  to  rarely  occur,  valid  results  may  be.,  expected* 
Areas  (3)  and  (4)  involve  the  use  of  unknown  parameters. 
Accordingly,  the  validity  of  results  from  investigating  these; 
two  areas  cannot  be  assured. 

The  search  for  mechanisms  peculiar  to  high  energy  sputtering 
was  made  using  the  (111)  surface.  This  surface  was  considered 
the  one  most  likely  to  show  additional  mechanisms  since  it  is 
the  most  densely  packed  surface.  The  failure  to  find  mechanisms 
(part  A  of  this  section)  peculiar  to  high  energy  sputtering 
was  not  too  surprising. 

The  sputtering  of  the  (110)  surface  by  ions  directed  Into 

(110)  channels  was  investigated  by  determining  the  sputtering 
ratio  for  ions  impacting  at  the  two  impact  points  located 
nearest  the  channels.  The  sputtering  ratio  curve  for  channel 
shots  Is  shown  In  figure  70.  The  discontinuity  between  the 
low  and  high  ion  energy  curves  is  assumed  to  be  the  result  of 
breakdown  of  the  potential  functions.  The  Important  result 
la  that  sputtering  always  occurred  for  channel  shots,  even 
when  40  keV  Ions  were  used. 

The  third  area  was  examined  using  xenon  to  sputter  the 

(111)  surface.  Argon-copper  potential  parameters  were  used 
aince  those  for  xenon-copper  were  not  available.  Runs  were 
made  at  3  and  5  keV.  At  both  energies,  the  ion  penetrated 

at  least  to  the  fourth  layer  of  atoms*  no  Indication  ofm  omen  turn 
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reversal  of  the  ion  was  observed.  The  Ion  was  found  to  be 
channeled,  just  before  penetrating  the  third  layer,  into  (110) 
channels.  It  initiated  cascades  which  propagated  into  the 
crystal. 

The  fourth  area  was  investigated  to  see  if  the  effect  of 
surface  geometry  on  the  sputtering  deposit  patterns  could  be 
determined  by  a  comparison  of  (111)  sputtering  of  a  face 
centered  cubic  crystal  with  sputtering  the  (0001)  surface  of 
a  hexagonal  close-packed  crystal.  Haslguti,  Hanada,  and 
Yamaguchl^^)  have  sputtered  zinc  with  8  keV  argon)  the 
deposit  pattern  showed  an  outlined,  equilateral  hexagon  with  • 
a  central  haze.  An  attempt  was  made  to  reproduce  this 
pattern  with  the  simulation  model  modified  for  zine-zinc 
interactions.  The  potential  parameter  was  adjusted  for  a 
Born-Mayer  type  potential  V(r)  ■  where  A  -  5?. 

(Z^Zg)^^  keV.  This  relationship  was  determined  by  Ander¬ 
sen  and  Sigmund^?).  The  argon-zinc  potential  parameters 
(assuming  a  Born-Mayer  type  potential)  were  approximated  by 
using  those  for  argon-copper  interactions.  Neither  the 
sputtering  ratio  nor  deposit  pattern  matched  the  reported 
results.  This  was  not  unexpected  since  the  sine  atoms  In 
the  crystal  must  be  represented  as  ellipsiods  (with  the  minor 
axis  in  the  basal  plane)  rather  than  represented  as  spheres. 

The  (100)  surface  was  sputtered  at  100  ev  to  investigate 
the  possibility  of  a  45  degree  pattern  rotation  for  low  ion 
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energy  sputtering.  The  number  of  atoms  which  were  sputtered 
was  extremely  small.  All  sputtered  atoms  had  energies  less 
than  10  ev.  It  was  found  that  the  atoms  sputtered  most 
frequently  from  the  regular  surface  were  n.n.n.’s  (figure  71)  • 
Bombarding  vacancy  and  stub  surface  conditions  also  resulted 
In  a  predominance  of  n.n.n.’s  being  sputtered  rather  than 
sputtering  of  the  atoms  forming  the  square  about  the  target 
atom. 

E.  Results  Common  to  All  Surfaces. 

Momentum  focusing  was  observed  to  occur  only  into  the 
crystal  for  the  three  copper  surfaces  studied.  This  was 
particularly  evident  from  atom  track  displays  of  (111)  surface 
sputtering.  If  an  atom  underwent  momentum  reversal,  it  was 
always  by  reflection  from  atoms  located  no  deeper  in  the 
crystal  than  the  third  layer.  Even  for  7  keV  ion  energies, 
crystal  atoms  located  deep  In  the  crystal  were  always  driven 
inward.  A  surprisingly  large  number  of  atoms  with  large 
energies  were  found  to  move  between  second  and  third  layers 
with  their  motion  nearly  parallel  to  the  surface  and  with  a 
small,  inward-directed  momentum  component.  The  energy 
which  appeared  to  be  delivered  in  impulses  to  atoms  still  at 
lattice  sites,  dissipated  through  the  crystal.  These  atoms 
rarely  acquired  more  than  20  ev  through  the  energy  impulse 
process.  Their  motion  was  restricted  since  they  were 
surrounded  by  other  atoms. 
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The  target  atom  in  the  (100)  and  the  (110)  surfaces, 
when  struck  at  near  zero  impact  parameter  by  the  ion, 
always  transmitted  the  majority  of  its  energy  billiard-ball 
fashion  in  close  packed  rows  perpendicular  to  the  surface. 
(Atoms  to  which  this  energy  was  transferred  escaped  through 
the  back  face  of  the  crystal.)  There  were  a  few  cases  In 
which  one  of  these  atoms  would  initiate  another  chain.  The 
chain  was  never  more  than  2-3  atoms  in  length  before  the 
energy  was  either  wholly  dissipated  or  the  chain  stopped  by 
divergence  of  the  momentum  to  form  numerous  small  cascades 
directed  into  the  crystal  interior . 

The  unique  arrangement  of  atoms  in  the  first  layer  was 
observed  to  be  a  dominant  factor  in  determining  ejection 
directions  for  all  surfaces.  Atoms  which  were  ejected  other 
than  nearly  normal  to  the  surface  were  always  influenced  by 
their  neighbors.  Atoms  which  were  ejected  at  an^es  near 
45  degrees  to  the  normal  would  reflect  from  their  neighbors, 

often  ending  up  as  normally  ejected  atoms.  This  effect  was 

.  « » 

most  pronounced  for  (111)  sputter'ng  since  the  six  atoms 

acted  as  a  lens,  but  it  was  also  seen  for  (110)  sputtering 
in  which  a  lens  is  formed  by  second  layer  atoms.  The  (100) 
surface  showed  a  strong  lens  effect  although  one  would  not 
necessarily  assume  that  it  occurs.  Sputtered  atoms, 
originally  adjacent  in  the  (100)  surface,  were  often  observed 
to  be  ejected  almost  simultaneously.  In  these  situations  they 
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were  ejected  nearly  normal  to  the  surface.  This  may  be  the 
genesis  of  dimers  recently  observed  by  Woodyard^-^. 

The  spot  regions  of  the  deposit  patterns  from  (100)  and 
(111)  sputtering  always  contained  more  high  energy  (^10  ev) 
than  low  energy  atoms.  But  the  energy  distribution  was  more 
uniform  than  increasing  with  the  ion  energy.  It  was  not 
peaked  at  any  one  or  group  of  energies.  Atoms  with  low  Ion 
energies  were  predominantly  found  in  the  central  region  of  the 


pattern. 


t 
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6.  Conclusions . 


The  observable  quantities  of  sputtering  appear  to  be 


Interdependent  on  only  the  sputtering  mechanism*.  There  Is 
no  evidence  that  these  quantities  are  interrelated.  The 
qualitative  and  quantitative  data  from  the  simulation  indicate 
that  the  deposit  pattern,  sputtering  ratio*  and  energy 
distribution  of  sputtered  atoms  cannot  be  correlated  with 
each  other;  one  cannot  predict  a  sputtering  ratio  from  an 
energy  distribution.  But,  each  of  these  quantities  can  be 
cross-correlated  between  surfaces. 

The  main  features  of  deposit  patterns  appear  to  be 
determined  only  by  the  surface  geometry  of  the  crystal.  The 
formation  of  (111)  pattern  spots  Is  attributed  primarily  to 
aaslated  focusing  by  the  hexagonal  lens.  The  predominance  of 
trigonal  rather  than  hexagonal  symmetry  is  considered  a 
natural  result  of  the  brief  (112)  channeling  observed  In  the 
mole  mechanism.  The  appearance  of  streaks  between  the 
spot  pairs  and  the  appearance  of  a  hexagonal  spot  pattern  at 
Z  keV  suggests  that  use  of  a  hemispherical  collector  in  the 
simulation  will  show  the  three  (Ilk)  epote  which  were  observed 
by  Southern  and  Kobinaon* 33  * .  However,  no  distinction  can 
be  made  between  the  (110)  and  (1H)  epota  on  the  basis  of 
pattern  location.  If  the  mole  mechanism  le  predominant  at 
certain  energies  the  (110)  spots  should  be  more  intense  than 
the  (Ilk)  epote.  If  the  souse* e  mechanism  le  predominant. 
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the  (110)  and  (114)  spots  should  have  more  equal  intensities. 

The  nearly  exact  simulation  of  a  (100)  pattern  at  3  keV 
is  considered  one  of  the  best  arguments  that  sputtering  is 
mostly  a  surface  phenomenon.  Neither  the  potential  form 
nor  its  parameters  for  copper-copper  and  argon-copper  inter¬ 
actions  are  known  with  certainty.  If  sputtering  was  a  deep 
phenomenon,  one  might  accept  an  argument  that  spots  in  the 
pattern  could  be  produced  without  exact  knowledge  of  the 
potential.  However,  the  hypocycloid  outline,  definitely  present 
in  the  simulation  pattern,  is  considered  the  feature  of  the 
pattern  which  would  be  most  sensitive  to  small  variations  in 
the  potential  if  sputtering  involved  more  than  the  first  few 
layers.  This  conclusion  is  consistent  with  the  apparent 
rotation  of  the  1  keV  pattern.  Atoms  found  in  ^211^  regions 
were  sputtered  when  the  energy  transferred  by  the  ion  to  the 
target  was  as  small  as  40  ev.  This  compares  favorably  with 
a  calculated  transfer  value  of  75  ev  for  a  0.5  A  impact  - 
parameter  a  sputtering  threshold  energy  of  50  ev 

has  been  reported  ^32)t 

The  inability  to  produce  consistent  ( 100 )  patterns  at  a 
given  binding  energy  for  5  and  7  keV  sputtering  is  undoubtedly 
due  to  the  unknown  condition  of  the  surface.  It  would  be 
unreasonable  to  assume  that  the  binding  energy  is  a  decreasing 
function  of  Ion  energy  only  for  the  (100)  surface.  The 
disruption  of  the  surface  by  the  first  group  of  incident  ions 
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is  a  factor  which  cannot  be  ignored.  The  sputtering  ratios 
determined  for  the  (100)  surface  as  well  as  for  the  ( 1 10 ) 
and  (111)  surfaces,  were  within  reasonable  limits  when  the 
assumed  binding  energy  is  in  the  range  1.50  to  3*50  ev.  No 
definite  value  of  binding  energy  can  be  determined  from  the 
patterns  and  sputtering  ratios  unless  a  weighting  factor  is 
used  for  the  surface  condition  at  the  time  of  ion  impact. 

The  assignment  of  weighting  factors  would  be,  at  best,  a 
guess . 

Three  surprising  features  were  observed  in  simulation 
deposit  patterns.  The  first,  the  apparent  rotation  of  the 
(100)  pattern  at  1  keV,  will  be  further  discussed.  .  The 
rotation  cf  the  (100)  pattern  for  low  energy  sputtering  has 
not  been  observed.  The  results  of  the  simulation  indicate 
that  such  a  rotation  Is  possible.  This  conclusion  is  based  on 
empirical  rather  than  theoretical  considerations.  The  hypocycloid 
outline  found  at  3  keV  In  the  simulation  appears  to  be  formed 
by  the  same  atoms  which  form  (£11)  spots  at  1  keV. 
Furthermore,  no  explanations  have  been  previously  proposed  as 
to  why  the  outline  is  a  hypocycloid  In  experimental  patterns. 

One  would  expect  that  a  circular  or  perhaps  a  non-distinct 
outline  would  be  observed  when  the  (100)  surface  is  sputtered. 
Accordingly,  It  is  proposed  that  the  hypocycloid  outline  is  a 
result  of  the  inability  of  the  (100)  surface  to  completely 
focus  atoms  sputtered  by  a  squeeze  mechanism  for  bombardment 
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energies  at  which  the  energy  propagation  is  not  confined  in  the 
surface. 

The  second  feature  was  the  unique  form  of  the  bonds  which 
connected  the  spots  in  the  simulation  pattern  for  (111) 
copper  sputtered  at  5  keV.  The  appearance  of  the  similar 
outline,  a  spherical  triangle,  was  observed  in  a  pattern  of  (ill) 
copper  sputtered  with  1.5  keV  krypton  by  Yurasova  and 
Bukhanov^  34) ,  The  third  feature  was  the  hexagon  spot  pattern 
at  2  keV,  a  pattern  seen  by  Anderson  and  Wehner^ for 
(111)  copper  sputtered  with  400  ev  mercury.  The  similarities 
observed  in  the  simulation  patterns  with  these  anomalous 
features  of  the  experimental  patterns  suggested  that  common 
factor  other  than  the  (111)  copper  surface  might  be  present. 
The  momentum  ratios  of  400  ev  mercury-2  keV  argon  and. 
L5keV  krypton-5  keV  argon  are  1.0  and  0.8  respectively. 

This  is  not  considered  to  definitely  establish  a  sputtering 
correlation  based  only  on  ion  momentum.  The  momentum  ratios 
for  1.5  keV  krypton  and  4  and  3  keV  argon  are  0.9  and  1.0 
respectively.  Ideally,  the  spherical  triangle  would  be  observed 
in  the  3  keV  (111)  simulation  pattern.  Additionally,  when 
(111)  copper  is  sputtered  at  high  temperatures,  the  triangular 
outline  has  been  observed  to  become  more  pronounced  ^34)  , 

The  absence  of  such  a  momentum  scaling  effect  could  be 
easily  shown  by  experimentally  sputtering  (111)  copper  with 
Ion  momentum  comparable  to  that  of  400  ev  mercury. 
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The  general  consistency  of  the  numerical  values  of 
sputtering  ratios,  and  energy  distributions  of  sputtered  atoms 
and  likeness  of  the  patterns  with  experimental  data  is  a 
strong  argument  for  the  validity  of  a  computer  simulation  of 
sputtering.  It  is  remarkable  that  the  consistency  is  as 
good  as  it  is.  The  model  uses  only  a  repulsive  potential,  the 
crystal  size  used  in  the  simulation  is  an  infinitesimal  portion 
of  the  smallest  laboratory  specimens,  and  the  potential  form 

4 

and  parameters  are  comparatively  crude.  It  is  concluded 
that  the  results  obtained  substantiate  the  concept  of 


transparency  and  the  occurrence  of  momentum  focusing  onl; 


within  the  surface  layer. 
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APPENDIX  A 


The  Beam  Model 

A.  Impact  areas  and  Impact  points. 

Each  crystal  surface  contains  an  intrinsic ,  plane  geometric 
shapej  hexagon  for  (ill)  surface,  square  for  (100)  surface, 
and  rectangle  for  (110)  surface.  A  volume  element  of  the 
(hkl)  surface  is  defined  by  the  area  of  this  intrinsic  shape 
and  a  depth  of  some  number  of  (hkl)  planes.  This  volume  is 
chosen  so  that  by  translation  along  axes  of  a  Cartesian 
coordinate  system  the  entire  crystal  may  be  generated.  The 
smallest  intrinsic  area  which  can  be  chosen  for  each  surface 
and  still  satisfy  the  translation  requirement  is  shown  in 
figure  72,  A  finite  number  of  points  Is  symmetrically  distributed 
within  each  of  these  areas  to  represent  the  infinite  set  of 
possible  points  of  Impact  for  an  incident  ion.  These  smallest 
areas  are  further  divided  into  representative  Impact  areas 
which  are  seen  from  figure  73  to  be  degenerate  under  appro¬ 
priate  rotations  and/or  inversions  of  the  coordinate  axes. 

Since  the  sectors  are  degenerate,  the  set  of  impact  points 
In  each  sector  is  also  degenerate!  only  the  points  contained 
in  one  Impact  area  need  be  used  to  represent  bombardment  of 
the  entire  area.  The  coordinate  axes  rotation  and  inversion 
schemes  are  discussed  in  Appendix  C  (Deposit  Pattern 
Production) . 

The  independence  of  the  Impact  point  set  used,  with 


respect  to  number  and  kinetic  energy  of  atoms  sputtered,  was 

tested  using  two  sets  of  points.  No  dependence  was  found. 

Results  for  the  (111)  surface  using  a  3  keV  ion  are  shown  in 

figure  74  as  an  example.  Set  1  points  are  those  shown  in 

figure  73,  set  2  (not  shown)  is  a  set  of  eleven  points  which 

are  located  between  the  points  of  set  1. 

B.  Neutralization  of  Beam  ions. 

The  assumption  that  argon  ions  are  neutralized  prior  to 

impact  on  copper  may  be  inferred  from  results  of  a  theoretical 

(27) 

study  of  secondary  electron  emission  by  Harrison  et.  al.  . 

Consideration  of  atom-atom  rather  than  ion-atom  interactions 

gave  results  in  reasonable  agreement  with  experimental  data. 

A  supporting  argument  for  neutralization  is  based  on 

Hagstrum's  theory  of  Auger  ejection  of  electrons' J  .  The 

probability  of  an  ion  being  neutralized  In  dx  at  x  Is: 

Ft  (x,v)  «a  exp  |  -  exp  [-a(x-xm)J  -a(x«xm)^ 

where  xm  ■*  (1/a)  In  (A/av)  is  the  value  of  x  where  Ft  is  a 

maximum.  The  parameters  A  and  a  occur  in  the  transition 

rate  function,  and  v  is  the  velocity  of  an  ion  starting  at  x  -  oo. 

Hagstrum  used  tungsten  as  an  example  and  obtained  the  value 

of  A  by  empirical  means,  the  value  of  a  from  published  data. 

Rather  than  attempt  an  exact  proof  of  neutralization  for  argon 

on  copper,  Hagstrum's  results  for  argon  on  tungsten  have  been 

used  to  give  an  order  of  magnitude,  at  worst,  approximation. 

o 

Accordingly,  l.o3^xm{2.18  A  for  1  to  10  keV  argon  Ions  on 
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copper |  this  is  a  reasonable  distance  from  the  crystal  surface 

•  » 

and  one  may  assume  that  neutralization  occurs. 

C.  Ion  Deflection  by  Surface  Potential. 

As  a  first  approximation,  it  is  assumed  that  a  singly 
positive-charged  argon  ion  sees  the  crystal  as  a  single,  fixed 
scattering  center  having  2-1  positive  charge.  The  well  known 
equation  for  a  central  force-induced  hyperbolic  orbit  is  used 
in  conjunction  with  figure  75  to  determine  the  deflection. 

r(0)  -  (L2/mk)  /  j^-1  +  (l  +  ZEL^mK2)1^  cos(  0-  Q*  )J 

L  “  ion  angular  momentum  m  «  ion  mass 
E  -  ion  energy  K  ■ 

Q  -  angle  of  closest  approach 
The  angular  momentum  L  is  determined  at  r-*oo,  C-+  W 
to  be  L  ■  mvb  where  m v^/2  **  E  and  b  is  the  impact  parameter. 
The  equation  for  r(  0  )  is  put  in  a  more  convenient  form  by 
the  substitution  L2  -  ( mvb )  *"  ■  2mEb2  to  give 
r  (©)  -  b  (2Eb/k)/£  -1  ♦  [l  +  (2Eb/k)2Jl/2  cos(0-ft,  ) 
The  angle  of  closest  approach  A,  is  determined  as  a  function 
of  the  product  Eb  by  requiring  the  denominator  to  vanish  for 
6*T  .  Once  ft  a  function  of  Eb)  is  known,  the 
deflection  ratio  r/b  is  calculated  at  the  crystal  surface. 

The  deflection  calculated  is  for  an  unneutralized  ion,  thus 
greater  than  that  for  an  ion  which  is  neutralised  at  some 
distance  in  front  of  the  surface.  Plots  of  r/b  at  the  surface 
and  ft-  functions  of  the  ion  kinetic  energy-impact  parameter 
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product  are  contained  in  figure  76.  The  numerical  values  are 
+  + 

for  an  Ar  — •  Cu  system. 

A  more  sophisticated  approximation,  scattering  from  a 

fixed  dipole  of  same  charge  sign,  requires  a  messy  Integration. 

Rather  than  follow  this  line,  one  may  expect  from  the  nature 

of  the  problem  that  the  path  of  the  approaching  ion  will  undergo 

some  oscillatory  motion  or  perhaps  corkscrew  motion  as  the  ion 

is  influenced  by  the  surface  potential.  In  either  case  one 

would  expect  that  the  net  acceleration  of  the  ion  parallel  to 

the  crystal  surface  would  be  no  greater  than  that  due  to  a 

single  fixed  scattering  center.  With  these  considerations, 

figure  76  is  used  to  determine  the  percentage  of  ions  which 

will  be  appreciably  deflected.  As  an  example  of  an  appreciable 

ion  deflection,  consider  r/b  -  2.0 j  from  figure  76,  Eb  ■  3.5  x 
-2  ° 

10  keV-A.  An  energy  range  of  1-10  keV  for  argon  ions, 

that  used  in  this  study,  corresponds  to  a  range  of  Impact 

o 

parameters  of  0.035  -  0.0035  A.  The  fraction  of  surface 
area,  and  therefore  fraction  of  ions  which  will  be  appreciably 
deflected  is  (0.035)2/( 1.26) 2  -  7.7  x  10"4  for  1  keV  ions  and 

c  0 

7.7  x  10”°  for  10  keV  Ions  (1.26  A  is  the  copper  atom 
effective  radius  in  the  crystal).  This  negligibly  small  fraction 
of  ions  cannot  Influence  the  macro  aspects  of  the  sputtering 
model. 

D.  Equilibrium  State  During  Impact. 

The  use  of  a  single  atom  approximation  to  the  beam  with 
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the  crystal  model  described  In  the  main  text  requires  that  the 
crystal  region  be  in  an  equilibrium  state  at  the  time  of  an 
ion  impact.  It  is  not  necessary  that  this  equilibrium  state 

t 

be  identical  to  previous  equilibrium  states  since  the  variety 
of  surface  configurations  available  in  the  model  provides  for 
random  surface  conditions.  It  is  necessary  that  the  time 
required  for  the  crystal  region  to  return  to  an  equilibrium 
state  be  small  with  respect  to  the  arrival-time  intervals  of 
the  lens.  Satisfaction  of  this  condition  is  determined  by 
comparing  an  experimental  beam  flux  over  the  area  of  the 
crystal  face  of  the  model  used  to  the  time  required  for 
completion  of  all  energetic  collisions  in  the  model. 

„2 


Beam  intensities  of  the  order  of  10^* amps/ cm^  were  used 


by  Magnus  on  and  Carlston 


(12) 


If  a  beam  of  this  intensity 
is  Incident  on  a  crystal  surface  area  of  less  than  lCpA^  such 


as  In  the  model,  the  Ion  flux  over  this  area  is  less  than  100 
ions/sec,  an  ion  arrival-time  interval  of  10“2  seconds)  all 
energetic  collisions  in  the  model  are  completed  within  about 
10  seconds.  The  relaxation  tiros  of  the  region  is  thus  much 
smaller  than  the  ion  arrival- time  intervals  and  the  condition, 
that  the  crystal  region  be  in  an  equilibrium  state  at  the  time 
of  impact,  is  satisfied. 
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APPENDIX  B 

Positioning  the  Ion  v. 

The  ion  is  positioned  tangent  to  the  first  target  atom  it 
will  strike.  Since  this  is  a  dynamic  rather  than  static  process 
it  is  not  necessary  that  this  be  a  stable  position  on  the  crystal 
surface.  Figure  77  shows  the  ion  at  its  arbitrary  initial 
position  and  calculated  final  position,  both  with  respect  to  the 
impact  point  and  target  atom.  The  initial  position  is  a  small 
distance  in  front  of  the  surface,  beyond  the  eroded  potential 
range  of  the  crystal  atoms.  This  position  is  described  by  a 
vector  originating  at  the  impact  point  and  having  direction 
parallel  but  opposite  to  the  ions  velocity  vector.  The  vector 
rg,  from  the  impact  point  to  the  target  atom,  is  known 
since  the  target  atom's  coordinates  are  known;  the  vector  r^ 
is  known,  and  it  is  desired  that  vector  r^g  have  a  magnitude 
equal  to  the  distance  between  centers  of  two  tangent  atoms, 
2rQ.  The  law  of  cosines  is  used  to  compute  the  magnitude  of 
r^  which  lies  along  r^.  Accordingly,  the  following  calculations 
are  made* 

cos  Ot  «  ?  .  r 

irii  TH 

(r~)k  -  (2 r  coset  )  (r  )  +  (r^  -  (2r  )3)  »  0 
“  “  3  2.0  — 

< 2-3 )  -  Tg  coa*4  +  jr|  cos2e<.  -  (r|  -  (2r0)2)] 

The  positive  square  root  solution  is  chosen  to  give  the  tangent 
on  tie  outside  hemisphere  of  the  target  atom. 
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The  ion  positioning  is  accomplished  automatically  for  each 
run  by  subroutine  START  t  which  is  contained  in  the  computer 
program  for  the  sputtering  simulation. 
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APPENDIX  C 

I 

Production  and  Analysis  of  Sputtering  Deposit  Patterns 
A.  Production. 

A  sputtering  deposit  pattern  represents  the  intersection 
points  of  atoms’  velocity  vectors  with  the  surface  of  a 
collector  plate.  The  sputtered  atoms  from  the  simulation 
are  collected  on  a  flat  plate  by  determining  these  points  of 
intersection.  Each  atom  which  exits  through  the  crystal 
surface  is  initially  assumed  to  have  been  sputtered  and  a  data 
card  has  been  prepared  for  each  one.  Data  of  particular 
interest  are  each  velocity  component  magnitude  and  the  kinetic 
energy  perpendicular  to  the  surface}  data  of  secondary 
Interest  are  atom  number,  impact  point  used  and  ion  kinetic 
energy.  (The  use  of  these  last  data  will  be  discussed  shortly.) 

It  is  recognized  that  each  sputtered  atom  has  lost  some 
energy  to  overcome  the  surface  binding  energy.  An  assumed 
value  of  binding  energy  is  subtracted  from  the  perpendicular 
kinetic  energy  and  a  new  perpendicular  velocity  component  is 
calculated.  Parallel  velocity  components  are  normalized  to  the 
new  perpendicular  component  to  give  a  two  dimensional  coordin¬ 
ate  point.  This  point  is  the  intersection  of  the  atom's 
velocity  vector  with  an  imaginary  collector  plate  placed  at  unit 
distance  from  the  target  surface.  Each  point  is  then  rotated 
and/or  mirrored  about  the  coordinate  axes  to  give  the  inter¬ 
section  of  a  velocity  vector  which  would  have  resulted  for  an 
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ion  impacting  at  the  corresponding  impact  point  In  each  of 
the  other  impact  areas.  The  impact  areas  are  shown  in 
figure  73  and  the  coordinate  point  rotation  and  mirroring  values 
are  listed  in  program  DATASORT  which  is  used  to  generate 
the  points.  Bach  point  is  plotted  using  program  DATAPLOT 
with  a  CDC160A  computer  and  a  CaiComp  plotter  to  give  a 
point  plot.  The  dimensions  of  the  plots  are  In  target-to- 
collector  or  T-O  units  since  these  points  have  been  normalised 
to  unit  target- to-eollector  distance.  The  scale  which  has 
been  used  permits  plotting  deposit  points  of  atoms  which  have 
been  sputtered  within  an  escape  cone  of  about  63  degrees 
(57  degree  cone  shown  in  figures).  This  has  been  found 
satisfactory  to  contain  all  pattern  features  of  interest. 

The  point  plots  show  only  point  patterns  and  therefore 
do  not  accurately  simulate  experimental  patterns  which  are 
area  density  patterns.  The  conversion  of  a  point  plot  to  a 
smooth  area  density  pattern  is  made  by  photographing  the  point 
plot  with  the  camera  defocus ed  such  that  no  single  point  is 
distinguishable  but  high  and  low  density  areas  are  prominent. 
Developing  and  printing  is  controlled  to  bring  out  the  high 
density  areas  while  maintaining  the  hase  background.  Loss  of 
intensity  in  some  spot  regions  cannot  be  avoided  such  as  seen 
in  figures  6b  and  12b.  The  process  must  be  adjusted  for 
each  pattern i  the  sequence  of  photography,  developing,  and 
the  printing  of  the  positive  image  is  highly  dependent  on  the 
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ability  of  the  photographer.  Concurrently,  the  quality  of 
reproduction  of  these  patterns  in  printed  form  depends  on  the 
plate  preparation,  the  printing  press,  and  the  paper  used.  s  < 

B.  Analysis,  ;  . 

A  30  x  30  square  grid  is.  placed  over  the  central  3,0  x  3*0 
unit  square  of  the  raw  pattern  by  program  DATAGRID. 

The  identification  number  of  an  atom  and  the  impact  point  and 
Impact  area  of  the  run  in  which  that  atom  was  sputtered  are 
recorded  for  each  atom  in  the  grid  square.  The  total  number 
of  atoms,  the  total  energy,  and  the  average  energy  per  atom 
for  each  grid  square  are  printed  in  separate  square  arrays. 
The  individual  grid  square  data  provides  for  correlation 
between  an  atom’s  crystal  location  and  its  deposit  point  in 
the  pattern)  the  square  arrays  of  the  number  of  atoms  and 
the  total  and  average  energy  densities  may  be  compared 
directly  with  the  point  plot  or  smooth  pattern  for  analysis  of 
pattern  features. 
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APPENDIX  D 


Section  1 


Section  2 


Section  3 


Memory  block  allocations  and  functions  are 
established.  All  storage  cells  are  zeroed. 

Input  and  output  formats  are  specified. 

Constants  are  set,  and  the  target  material, 
ion  species,  and  crystal  face  to  be  sputtered 
are  read  in. 

The  ion  energy,  impact  area,  and  impact 
point  are  read  in.  Constants  peculiar  to  the 
run  are  established.  The  appropriate  sub¬ 
routine  is  called  to  build  the  crystal,  and 
crystal  boundaries  are  calculated.  Subroutine 
START  is  called  to  position  the  ion.  The 
initial  value  of  At  is  calculated.  Initial 
coordinates  of  all  atoms  are  assigned,  and 
their  velocity  components  are  zeroed  {except 
for  the  ion).  The  initial  coordinates  of  the 
crystal  atoms  are  printed. 

Forces  are  calculated  by  calling  subroutine 
STEP.  Atoms  are  moved  to  their  intermediate 
positions  in  the  two  step  cycle.  Subroutine 
STEP  is  called  again,  and  atoms  are  moved 
to  their  final  positions.  Final  velocities  are 
then  computed.  Force  components  are 
zeroed  in  preparation  for  the  next  tlmeetep. 


The  maximum  kinetic  energy  ie  determined  for 
the  calculation  of&t. 

Section  4s  The  time  remaining  until  cutoff  ie  determined. 

If  there  ia  insufficient  time  to  complete 
another  tlmestep,  terminal  data  ie  printed. 
Potential  energy  is  calculated  and  summed 
with  kinetic  energy  to  give  the  total  energy 
for  energy  balance  check  ( manual  check) .  The 
data  for  atoms  having  potential  energy  greater 
than  the  thermal  energy  are  printed. 

Section  5  s  Atoms  which  have  kinetic  energy  but  not 

potential  energy  are  assumed  to  be  free  of 
the  crystal.  They  are  assigned  L»CUT*»1. 

The  surface  through  which  an  atom  exited  le 
determined  and  a  code  assigned.  A  maximum 
potential  energy  is  found  (for  ht  calculation) 
among  atoms  which  do  not  have  DCUT-l.  If 
the  maximum  potential  energy  ie  less  than  a 
minimum  value  the  terminating  process  begins. 

If  this  energy  is  greater  than  the  minimum 
value*  another  tlmestep  begins. 

Section  6s  Pertinent  data  for  all  atoms  ie  printed.  Atoms 
which  have  exited  through  the  front  of  the 
crystal  or  will  exit  through  the  front  are 
assumed  to  bs  sputtered.  Data  for  theee 
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atoms  areprinted  and  a  data  card  punched 
for  each  atom.  A  new  data  card  is  then 
read  Into  initiate  another  run. 
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APPENDIX  E 
Glossary  for  FCCSPUT 

AC  Parameter  for  target  force  function  ccrr action. 

AIX 

AIY  Floating  point  form  of  IX ,  IY,  12, 

AI2 

ALFA  Cosine  of  the  angle  between  vectors  Rl*  R2. 

ALFA2  ALFA  squared. 

B  Reciprocal  of  magnitude  of  atom  velocity. 

SENGY  Energy  which  an  atom  within  the  crystal  at  shut¬ 

down  must  have  to  be  considered  sputtered. 

BULLET  Variable  representing  primary  material. 

BX  Unsealed  x,  z  coordinates  of  the  impact  point 

BZ 

* 

CELS  Frictional  force  multiplier.  (See  CVS) 

COX 

COY  Direction  cosines  of  primary  velocity  vector. 

CO  2 

COXI 

COYI  Negative  values  of  COX,  COY,.  CD2. 

C02I 

CVB  A  constant. 

CVD  Converts  meters  to  angstrom  units, 

CVE  Converts  electron  volts  to  joules. 

CVM  Converts  atomic  mass  units  to  kilograms. 

CVR  Converts  lattice  units  to  angstrom  units. 

DCOX 

DCOY  Direction  cosines  of  sputtered  atom  velocity  vector. 

DC02 

DFF  Distance  difference  between  nearest  neighbor 

distance  and  actual  atom  ^differences. 

DIST  Distance  between  any  two  atoms. 
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DRX 

DRY  x,  y,  z  components  of  DIST. 

DR3 

DT  Length  of  timestep  to.  aeconda. 

DTI  Number  of  lattice  units  most  energetic  atom 

may  move  in  one  timestep. 

DTOD  DT/CVD  -»  a  ratio  used  to  avoid  repeated  division. 

DTOM  DT/PTMAS  -  a  ratio  used  to  avoid  repeated 

division. 

DTOMB  DT/PGMAS  -  a  ratio  used  to  avoid  repeated 

division. 

DX 

DY  x,  y,  z  distances  atom  has  moved  from  initial 

DZ  position. 

ECUT  A  lower  limit  on  an  atom's  potential  energy.  If 

energy  is  less  than  or  equal  to  ECUT  the  program 
shuts  down. 

El  A  cutoff  energy. 

EMAX  The  maximum  energy  encountered  to  any  cycle. 

EV  Primary  energy  In  electron  volts. 

EVR  Primary  energy  in  kilo-electron  volts. 

EXA  Potential  function  parameters. 

EXB 

FA  The  component  force  increment  on  an  atom. 

FAC  The  minimum  distance  the  primary  is  positioned 

in  front  of  the  first  xz  plane  at  start  time. 

FM  A  small  number  used  In  checking  potential  energy 

zero  point. 

FOD  FORCE/DIST  -  a  ratio  used  to  avoid  repeated 

division. 

FORCE  Numerical  value  of  the  force  function  with  a 

variable  parameter. 


FORF 

Target  atom  force  function. 

FFTb 

The  corrective  force  value  at  ROE, 

FPTF 

The  corrective  force  function. 

FRC 

Numerical  value  of  the  target  force  function  at 
ROE. 

FX 

FY 

x ,  y,  z  components  of  totai  force  on  an  etpm* 

FZ 

FXA 

Force  function  parameter. 

GMAS 

Target  atom  mass  (in  a.m.u. } 

HDTOD 

1/2  DTOD  -  a  ratio  used  to  avoid  repeated 
division. 

HDTOM 

1/2  DTOM  -  a  ratio  used  to  avoid  repeated 
division. 

HDTOMB 

1/2  DTOMB  -  a  ratio  used  to  avoid  repeated 
division. 

HGMAS 

1/2  GMAS  -  a  ratio  used  to  avoid  repeated 
division. 

HTMAS 

1/2  TMAS  -  a  ratio  used  to  avoid  repeated 
division. 

ICUT 

Used  to  provide  output  prior  to  time  limit 
shutdown. 

IDUM 

Dummy  variable. 

IHB 

IHS  . 

IHT 

IH1 

Alpha-numeric  arrays  for  titling. 

IN 

» 

Odd-even  integer  used  to  determine  atom  site 
establishment. 

INDEX 

Integer  (0  or  1)  used  in  determining  dynamics 
cycle  step. 

I  NOW 

Time  program  has  been  running  in  seconds. 

IP 

ISHUT 

ISPX 

ISPZ 

ISS 

IT 

ITT 

IX 

IY 

12 

IXP 

IYP 

I2P 

JSHUT 

JT 

JTS 

JTT 

KCUT 

KT 

LCUT 

LL 

LS 

LSS 

M 


Subscript  value  of  atom.  Used  la  subi^tijje 
STEP. 

Time  left  prior  to  time  limit. 

Fixed  point  values  of  SPX,  SP2. 

Subscript  value  of  most  energetic  atom. 

Unsealed  fixed  point  x  coordinate  used  in  lattice 
generation.  Also  a  dummy  variable  in  function 
TIMEF. 

Odd-even  integer  used  to  determine  atom  site 
establishment. 

Number  of  x,  y,  z  planes  of  crystal. 

Crystal  dimensions  in  x,  y,  z. 

Cutoff  variable  based  on  total  potential  energy 
of  crystal. 

Unsealed  y  coordinate  used  in  crystal  generation. 
Variables  used  to  establish  atom  sites. 

Identifies  exit  point  of  atom. 

Unsealed  z  coordinate  used  to  establish  atom  site. 

Used  to  identify  acorns  which  are  not  included  in 
calculations. 

The  highest  numbered  atom  in  the  crystal. 

Sum  of  the  Miller  index  Integers. 

Used  to  identify  type  of  surface,  i.e.,  regular, 
stub,  vacancy. 

An  integer  used  to  begin  atom  numbering. 
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x 


N 

ND 

NPAGE 

NS 

N8HUT 

NT 

NTT 

PAC 

PEXA 

PEXB 

PFIV 

PFORF 

PFPTC 

PFPTF 

PFRC 

PFXA 

PGMAS 

PKE 

PKEY 

PDA 

PLANE 

PNUM 

POT 

POTF 

PPE 


Subscript  of  the  atom  to  be  removed  for 
Vacancy  surface* 

Data  output  Increment* 

Page  numbering  variable. 

Initial  print  statement  cycle  . 

Cutoff  variable  based  on  too  long  a  timeatep. 

Timestep. 

Times  tap  limit  before  shutdown. 

Same  as  AC  except  applicable  to  primary. 
Primary  force  function  parameters. 

A  constant  ■  0.5. 

Primary  force  function. 

PFPTF  evaluated  at  ROE. 

Primary  corrective  force  function. 

PFORF  evaluated  at  ROE. 

Primary  force  function  parameter. 

Primary  mass  in  kilograms. 

Kinetic  energy  of  an  atom. 

Y  component  of  kinetic  energy  of  an  atom. 
Crystal  plane  ( alphanumeric  variable). 

Same  as  PLA. 

Impact  point  (alphanumeric  variable). 

Potential  energy  between  two  atoms. 

Target  potential  function. 

Potential  energy  of  an  atom. 
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PPOTF 

PFTC 

PRI 

PTC 

PTE 

PTMAS 

QM 

QUIT 

RBX 

RB2 

ROE 

ROEM 


ROE2 


RX 

RY 

RZ 

RXBND 

RYBND 

R2BND 

RXI 

RYI 

RZI 

RXK 

RYK 

R2K 

RX8 

RYS 

RZ8 

RI 


Primary  potential  function. 

PPOTF  evaluated  at  ROE. 

Chemical  symbol  for  primary  material. 

POTF  evaluated  at  ROE. 

Totil  energy  of  an  atom  (potential  +  kinetic). 
Target  mass  in  kilograms. 

A  small  number  used  in  checking  kinetic  energy 
zero  point. 

Cutoff  variable  checked  against  total  potential 
energy. 

-Unsealed  x,  z  coordinates  of  impact  area 
reference  point 

Nearest  neighbor  distance. 

ROE  -  DTI  (one  timestep  distance  less  than 
n.n.  distance). 

ROE  squared. 


x,  y,  z  coordinates  of  atom  at  any  time. 


x,  y,  z  coordinates  of  crystal  boundaries  other 
than  zero . 

x,  y,  z  coordinates  of  atom’s  Initial  position. 

x,  y,  z  coordinates  of  temporary  position  of 
atom  during  force  cycle. 

x,  y,  z  coordinates  of  impact  point. 

Vector  from  Impact  point  to  initial  primary  position 
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R1X 

R1Y 

R1Z 

R1R2 

R2 

R2SQ 

R3 

SAVE 

sex 

SCY 

sez 

SLOW 

SPX 

SPZ 

88CZ 

SUR 

TAR 

TARGET 

TE 

TEMP 

THERM 

TI 

TIME 

TIMO 

TMAS 


x,  y,  z  coordinates  of  initial  primary  position* 

Scalar  product  of  vectors  Rl,  R2. 

Magnitude  of  vector  from  impact  point  to  first 
atom  hit  by  primary. 

R2  squared. 

Magnitude  of  vector  from  impact  point  of 
primary  start  position. 

1/2  POT. 

x,  y,  z  coordinate  scale  factors. 

Cutoff  variable  checked  against  a  long  DT . 

«  ■ 

x,  z  distance  from  Impact  area  reference  point 
to  impact  point. 

A  z  scale  factor  used  for  (ill)  plane  lattice 
generation. 

Plane  (alphanumeric  variable). 

Chemical  symbol  for  target  material. 

Target  material  (alphanumeric  variable). 

Total  energy  of  crystal  atoms  (kinetic  ♦ 
potential ) . 

Temperature  of  lattice  in  degrees  Kelvin. 
Thermal  energy  of  atom. 

Computer  time  program  has  been  running. 
Elapsed  problem  time. 

A  function  to  convert  seconds  to  minutes. 
Target  atom  mass  in  kilograms. 
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TFKE 


Total  kinetic  energy  of  crystal  atoms. 

*  *  y  x 

- t 

TPOT  Total  potential  energy  of  crystal  atoms. 

VOL  Magnitude  of  primary  velocity  vector. 

VSS  Storage  variable  for  velocity  components. 

VX 

VY  x,  y,  z  components  of  atom  velocity. 

V2 

X  Unsealed  x  coordinate  used  in  crystal  generation. 

XIT  Floating  point  form  of  IT  used  In  function  ,TIMO. 

Y  Unsealed  y  coordinate  used  in  crystal  generation. 

Z  Unsealed  z  coordinate  used  in  crystal  generation. 

2E  A  constant  **  0.0 

ZP 


Floating  point  form  of  JTT. 


Pattern  Production  and  Analysis  Frograrrs  (DATASORT 

DAT  APLOT,  DATAGPID) 
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Sputtering  frequency-location  dtagra 
for  (lOO)  surface. 
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(100)  Surface 
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Sputtering  frequency  -location  diagram 
for  ( 100)  surface. 
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Sputtering  frequency -location  diagram 
for  ( 100)  surface. 
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Sputtering  profits  of  atoms  sputtered 
from  (100)  surface. 


Sputtering  Ratio  Cotoms/mn) 


( 

I 


Ion  Energy  (Kev) 

Sputtering  Ratio  of  ( 100)  Surface 
Regular  Surface 

Figure  51a 
260 


Sputtering  Ratio  (atovns/ic 


4  6  6 

ion  Energy  (Kov) 

Sputtering  Ratio  of  (100)  Surface 
Vacancy  Surface 

Figure  51b 
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Spcttertng  frequency -location  dug  rare 
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(HO)  Surface 


Sputtering  frequency-location  diagram 
for  (110)  surface. 
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Figure  57 
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Sputtering  frequency-location  diagram 
for  (110}  surface. 
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Figure  58 
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(!!0)  Surface 


Sputtering  frequency-location  diagram 
for  (110)  surface. 
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Figure  61a 
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Figure  61b 
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Figure  62a 
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Figure  62b 
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Figure  63a 
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Figure  53b 
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Figure  64a 
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Figure  64b 
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sputtering  process  has  been  investigated  by  simulating;  the 
sputtering  of  single-crystal  copper  with  i-7  keV  ergon.  A  digital 
computer  was  used  to  build  the  crystal,  bombard  It,  and  move 
crystal  atoms.  Pour  mechanism#  observed  which  cause  surface 

atoms  to  sputter.  An  atom  is  sputtered  xvhen  (1)  it  is  squeezed 
out  of  ths  surface,  (2)  it  l#  scooped  out  when  another  atom  strikes 
Its  inner  hemisphere,  C3)  it  is  ejected  when  an  atom  passes  behind  it, 
and  (4)  it  is  knocked  out  fey  a  second  layer  atom  which  is  moving 
outward-  Nearly  all  sputtered  atoms  were  surface  atoms.  Second 
and  third  layer  atoms  were  sputtered  only  for  Ion  energies  greater 
than  5  keV.  They  were  sputtered  by  mechanisms  similar  to  the 
surface  atom  mechanisms,  ’’fifcllstoee  chains”  were  observed  to  fee 
directed  into  the  crystal,  and  momentum  focusing  was  observed  to 
cause  sputtering  only  When  it  occurred  in  close  packed,  surface  rows. 
Outward  directed  chains  were  not  observed.  Sputtering  deposit  '  " 
patterns,  sputtering  ratios,  and  sputtered  atom  energy  distributions 
were  obtained  for  (10Q),  (110),  and  (ill)  surfaces.  All  data 
compared  favorably  with  experimental  data. 
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